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Preface 

In the past years there has been considerable enthusiasm for the idea that nature might be 
super symmetric. The reasons are manyfold. Some people were motivated by the idea to unify 
bosons and fermions, others were attracted by the fact that local supersymmetry involves 
gravity. One of the most appealing features of supersymmetry (SUSY) is, however, that it 
provides a solution to the hierarchy problem by protecting the electroweak scale from large 
radiative corrections — provided that SUSY particles exist at or below the TeV scale. 

The experimental search for supersymmetric particles, which is one of the primary tasks at 
present and future colliders, relies on detailed phenomenological studies. In this thesis, I study 
the phenomenology of stops and sbottoms (the superpartners of the top and bottom quarks) 
in the framework of the minimal supersymmetric extension of the Standard Model (MSSM). 

This thesis is organized as follows: 

The first chapter serves as an introduction to the MSSM. I briefly review the model and explain 
why stops and sbottoms provide an interesting playground for exploring and testing this model. 

In the second chapter I discuss stop and sbottom production at lepton colliders. I give the 
analytical formulae for the production cross sections at e + e~ and ji + n~ colliders including 
(general) beam polarization. A numerical analysis shows the cross sections that can be expected 
at LEP2, an e + e~ Linear Collider, and a Muon Collider. In particular, I address the topics of 
what can be gained by using polarized beams and how one can determine the MSSM parameters 
by cross section measurements. 

The third chapter is devoted to stop and sbottom decays including 1-loop SUSY-QCD cor- 
rections. The calculations are done in the on-shell renormalization scheme using dimensional 
reduction. I present the analytical formulae for the 0(a s ) corrected decay widths and explain 
the subtleties that have to be taken into account in the calculations. Moreover, I perform a 
detailed numerical analysis of squark decay widths and branching ratios at tree level and 0(a s ). 
It turns out that stops and sbottoms can have very complex decay patterns. SUSY-QCD cor- 
rections can change the individual decay widths and branching ratios by a few ten percent. 
Hence they are important for precision measurements. 

This work grew from a series of workshops on physics at present and future colliders. Its 
aim is to give a comprehensive overview of the phenomenology of stops and sbottoms and to 
serve as a basis for Monte Carlo studies. 



ii 



Acknowledgements 

Above all, I thank my 'Doktorvater' Prof. Walter Majerotto for his guidance and advice, and 
the inspiring working atmosphere. He also gave me the opportunity to participate in several 
international workshops, including numerous visits to CERN. 

I thank Dr. H. Eberl for a close and lively collaboration. Likewise, I thank Prof. A. Bartl, 
Dr. T. Gajdosik, and Dr. W. Porod for many useful discussions. During the past two years I 
also enjoyed collaborating with Prof. K. Hidaka and Prof. Y. Yamada. I am grateful for this, 
too. 

I am indebted to the 'Rechengruppe' of our institute, led by Dr. G. Walzel, for excellent 
computing support. Moreover, I thank the Vienna CMS group for their hospitality during my 
visits to CERN. 

In a work that expands over years, dissappointments are of course not unknown. Though, joy 
clearly prevailed. In this context I wish to thank DI W. Adam, Dr. A. Boresch, DI B. Godel, 
Dr. R. Veenhof, and the members of the SKI Dojo Wien for their company and friendship. 

This work was supported in part by the "Fonds zur Forderung der wissenschaftlichen Forschung" 
of Austria, project no. P10843-PHY. 



Contents 

Preface 

Acknowledgements 

Contents 

1 The Minimal Supersymmetric Standard Model 

1.1 Definition of the Model 

1.2 Renormalization Gronp Equations 

1.3 Models of SUSY Breikmg" 

1.4 Electroweak Symmetry Breaking 

1.5 Sparticle Mixing and Mass Eigenstates 

|1.5.1 S quarks 

1.5.2 Sleptong 

1.5.3 Charginos 

1.5.4 Neutralinoi 

1.6 Interaction Lagrangiari 



1.6.1 Quark - Quark - Gauge Boson 








1.6.2 Squark - Squark - Gauge Boson 




1.6.3 Quark - Quark - Higgs Boson 








1.6.4 Squark - Squark - Higgs Bosor 


l 


1.6.5 Quark - Squark - Chargino . 




1.6.6 Quark - Squark - Neutralino 



iii 



CONTENTS 

1.6.7 Interactions with Gluinos| 

1.6.8 Squark - Squark - Gauge Boson - Gauge Boson 

1.6.9 Four-Squark Interaction 



2 Squark Production at Lepton Colliders 

271 Gross Sections for e + e~ Golliders . . . 



EXT LEP721 



2.1.2 Linear Golliders with Vs = 0.5-2 TeV 

2.1.3 Determination of Soft-Breaking Parameters — A Gase Study 
2.2 Gross Sections for Muon Golliders] 

■i SUSV-QCD Corrections to Squark Decays 

3~T Introduction 



3.2 General Aspects of SUSY-QCD Gorrections 
3.2.1 Renormalization 



3.2.2 Regularization: DREG vs. DRED 

3.2.3 Notations, Conventions, etc. . . . 
3.3 Self Energies to 0(aT> 

3.3.1 Squark Self Energy 

3.3.2 Quark Self Energy 

374 Subtleties 



|3. 4.1 Renormalization of the Squark Mixing Angle 
3.4.2 Renormalization of Mf, 



p.4.3 Renormalization of m a A~ a . . . 
3.5 Decays into Gharginos and Neutralinos 

p. 5.1 Vertex Correction! 

3.5.2 Wave-Function Correction . . . 



3.5.3 Renormalization of the Bare Couplings 



CONTENTS 

p. 5.4 Real Gluon Emission 
3.5.5 Numerical Results 



3TB Decays into W and Z Bosons 



3.7 Decays into Higgs Bosons 

3.8 Branching Ratios 

Open Questions 



Bibliography 



Chapter 1 



The Minimal Supersymmetric 
Standard Model 



1.1 Definition of the Model 

The simplest — and most popular — supersymmetric model is the straightforward supersym- 
metrization of the Standard Model (SM), where one introduces only those couplings and fields 
that are necessary for consistency. This is known as the Minimal Supersymmetric Standard 
Model (MSSM) |, |, The MSSM is build up as follows: 

• In addition to the gauge boson fields, spin-| "gaugino" fields are introduced. The 
partners of and are denoted B and W l . In analogy to the photon, the Z and 
the W ± bosons, one can form a photino 7, a Z-ino Z, and W /± -inos from the B and W l 
fields. The superpartners of the gluons are the gluinos g. 

• Quarks and leptons get spin-0 partners called SQUARKS and sleptons. As there has 
to be a superpartner for each degree of freedom, two bosonic fields are needed per SM 
fermion. They are named "left" and "right" states q L , q R and £ L , l R . 

• Moreover, one needs TWO complex Higgs doubletts with hypercharges ±1 in order 
to give masses to up- and down-type quarks and leptons, and to cancel anomalies. The 
Higgs fields are also assigned spin-| partners, the so-called higgsinos. 

The field content of the MSSM is shown in Table Supersymmetry (SUSY) [] in its local 
version includes gravity; the resulting theory is know as super gravity. The model then also 
includes the graviton (spin-2) and its fermionic partner the gravitino (spin-|). 

^or an introduction to supersymmetric field theory in general and the MSSM in particular see also [|], 0|; 
for a textbook on supersymmetry I refer the reader to B . 
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Table 1.1: Field content of the MSSM. 



Let us now turn to the Lagrangian of the MSSM. Clearly, gauge interactions are determined 
by the gauge group, which is the same as in the Standard Model: SU(3)c x SU(2) L x U(l)y. 
For a concise derivation of these interactions, see e.g. |2|, [|, [|, |f. Masses and couplings of the 
matter fields are determined by the superpotential W. The choice of the gauge group constrains 
W but does not fix it completely. Holding to the priciple of minimality, that means introducing 
only those terms that are necessary to build a consistent model, we get: 

W = E [(Me RihEj + (hnh H^D) + {h v ) l3 Q^U*} + fiH,H 2 (1.1) 

where i and j are generation indices. Contractions over SU(2) and SU(3) indices are understood. 
In particular, 

H X H 2 = HiaH%p = H^H® - HiH± (1.2) 

with e a/3 a totally anti-symmetric tensor used to contract over the SU(2)^ weak isospin indices 
a, (3 — 1, 2. Likewise, H X QD C = e afl ] H la Q a ^D c a where a = 1, 2, 3 is a colour index, etc.. The 3x3 
matrices hp, ha, and He are dimensionless Yukawa couplings giving rise to quark and lepton 
masses. Moreover, hp and hjj account for the mixing between the quark current eigenstates as 
described by the CKM matrix f5|. Notice that the same superpotential is obtained by requiring 
that baryon and lepton numbers be conserved (which is automatically fullfilled in the SM but 
not in the MSSM). 
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The Lagrangian derived from ( |1 . 1| ) is 

d 2 W 



SUSY 



E 



ipjipk + h.c. 



E 



;i.s) 



where fa are scalar and ipi fermion fields; W only depends on the scalar fields. The first term 
in ( |1.3|) describes masses and Yukawa interactions of fermions, while the second term describes 
scalar mass terms and scalar interactions. 

The interactions obtained in this way respect a symmetry called R-parity under which the 
"ordinary" fields (matter fermions, Higgs and gauge bosons) are even while their superpartners 
(sfermions, higgsinos and gauginos) are odd. As a consequence, all interactions involve an even 
number of SUSY particles ("sparticles"). This means that sparticles can only be produced in 
pairs, and any sparticle decay must lead to an odd number of sparticles. Hence in the MSSM 
the lightest supersymmetric particle (LSP) is stable. This leads to another important feature 
of the MSSM: Since the LSPs cannot decay some of them must have survived from the Big 
Bang era. Searches for so-called "exotic isotopes" have led to very stringent bounds which 
exclude any strongly or electromagnetically interacting stable particles with masses below a 
few TeV. The LSP of the MSSM must therefore be electrically and colour neutral ||. In turn, 
a neutral stable LSP is a good candidate for dark matter. For collider experiments this means 
that any decay chain of a sparticle will end in an arbitrary number of SM particles plus at 
least one LSP which escapes the detector, earring away some energy and momentum. In the 
context of the MSSM the typical SUSY signature is thus distinguished by missing (transverse) 
energy /momentum. Depending on different variants of the model, the LSP can be the lightest 
neutralino or the gravitino; the sneutrino has already been ruled out ||. 

We want to stress that R-parity is a symmetry that is somehow built in by hand due to our 
assumption of strict minimality. In principle, the following R-parity violating terms are allowed: 



XLLE C + \'LQD C + X"D C D C U C + [iH\L. 



;i.4) 



Here also generation indices have been suppressed. Within the MSSM breaking R-parity there- 
fore means breaking baryon and/or lepton number. This leads to a significantly different phe- 
nomenology compared to the R-parity conserving case. In particular, single sparticle produc- 
tion is possible, and the LSP is no longer stable. However, there are very stringent constraints 
on the couplings A, A', A", and //. If e.g., both baryon and lepton number were broken, this 
would lead to rapid proton decay. Together with other constraints |L0| and the wish to embed 
the MSSM in a Grand Unified Theory this requires that the couplings in ( |1.4| ) be at least very 
small, if not zero. However, in the way the MSSM is understood in this thesis, R-parity is 
conserved. (For the phenomenology of R-parity violating models see e.g. flQfl .) 

The scalar potential V is obtained from F— and D-terms: 



V = F*F i + \ 



D a D a + Jfjji + 



;i.s) 
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with 



1 - d-o) 



and 

D a = g 3 <P*T a <l> 1 D i = \g 2 <l>*(T i <f>, D' = \ gi Y^*<j> t . (1.7) 

T a = A a /2 (a = 1, ... 8) with A a the Gell-Mann-Low matrices. <t* (i = 1, . . . 3) are the Pauli 
matrices and Y^, = 2(Q — I 3 ) is the hypercharge. 

The Lagrangian as given in ( |1.3| ) conserves supersymmetry. However, in a realistic model 
SUSY must be broken Otherwise the masses of ordinary particles and their superpartners 
would be equal, which is not the case as we know from experiment. As the genuine mechanism 
of (dynamical) SUSY breaking is not yet understood, we parametrize it by inserting SUSY 
breaking terms by hand into the Lagrangian. The terms that break SUSY softly, i.e. do not 



induce quadratic divergencies, are fll 



gaugino mass terms — |M a A a A a , where a is the group index; 

scalar mass terms — M?.|0j| 2 ; 

trilinear scalar interactions A^ (pi(j)j(j)k] and 

bilinear terms —B^ (p^j + h.c. f\ 



They lead to the following explicit form of the soft SUSY breaking Lagrangian: 

-£ soft = \M X BB + \M 2 WW + ±M 3 gg + m 2 Hl {H^ 2 + m 2 H2 \H 2 \ 2 

+ M%\q L \ 2 + Ml\u c R \ 2 + Ml\d c R \ 2 + M\\i L \ 2 + Ml\e R \ 2 (1.8) 
+ (h E A E H x I L e R + h D A D H 1 q L d c R + h u A u H 2 q L u c R + B^H X H 2 + h.c.) . 

M x , M 2 , and M 3 are the U(l), SU(2), and SU(3) gaugino masses, respectively. m 2 H , m 2 H2 , and 
B\x are mass terms for the Higgs fields. The scalar mass terms Ml, M?, M?, M\, and M| 
are in general hermitean 3x3 matrices in generation space, while hjjAu, hnA^, and h E A E are 
general 3x3 matrices. Allowing all the parameters in ( |1.8|) to be complex, we end up with 
124 masses, phases and mixing angles as free parameters of the model. Notice that £ so ft also 
respects R-parity. Indeed, a R-parity violating term in Eq. ( |1.8| ) would lead to an unstable 
vacuum unless the same term also appears in the superpotential Eq. ( |1 . 1|) . 

2 In general also linear terms —C^i are allowed, where </>j is a gauge singlet field. 
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1.2 Renormalization Group Equations 



The physical quantities at the electroweak scale are related to their values at some high energy 
scale by renormalization group (RG) equations. We will not perform a RG analysis here but 
just add some quantitative arguments: 

One of the very appealing features of the MSSM is that it allows for gauge coupling unification 
at M x ~ 10 16 GeV jlj. The MSSM is thus compatible with a Grand Unified Theory (GUT). 
As in GUT models the gauginos all live in the same representation of the unified gauge group, 
gaugino masses are also unified at scales Q > M X - The 1-loop RG equations for the gauge 
couplings and the gaugino masses are 

with t = \a.{Q/Mx) and b a = 33/5, 1, —3 for a = 1, 2, 3, respectively |. One therefore has 



.I/, M 2 M 3 
9i 92 9s 

at any RG scale, up to small 2-loop effects. 



;i.io) 



We next consider the evolution of scalar masses. For simplicitly we here assume that the 
soft masses of squarks and sleptons are flavour-diagonal e.g., M| = diag(M? i , M?^, M?^). 
Moreover, we neglect Yukawa and trilinear couplings of the first and second generation, i. e. hjj = 
diag(0, 0, h t ), A v = diag(0, 0, A t ), etc.. The 1-loop RG equations for the masses of squarks and 
sleptons of the first two generations are then given by 

^ 2 jMl = -j:% g lCt\Mf (1.11) 

a=l 

with the sum running over the gauge groups. The are Casimir operators: Cf = |(-y) 2 for 
each scalar <p with hypercharge (Yq = |, Y u = — |, Y$ = |, etc.); C2 =3/4 (0) for <p = Q,L 
(U c , D c , E c ); and Cf = 4/3 (0) for <f> = Q, U c , D c (L, E c ). The right side of Eq. ( jETT] ) is strictly 
negative, so M? grows when being evolved down to the low scale. Moreover, owing to SU(3) 
contributions, squark masses grow faster than slepton masses. 

The mass-squared parameters of the third generation squarks and sleptons as well as of the 



Higgs fields also obey fll.ll ) but get additional contributions from Yukawa and trilinear cou- 



3 The coefficients b a in the MSSM are different from those in the SM, where b a = {41/10, -19/6, -7}, due 
to the richer particle spectrum in the loops. 
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plings: 



and 



16vr 2 - M? 3 = X, + X, - - g\ |M 3 | 2 - 6gj \M 2 \ 2 - - g\ \M,\\ (1.12) 

167T 2 | M? s = 2X, - |M 3 | 2 - | , 2 \M,\\ (1.13) 

16tt 2 | M? 3 = 2X b -^gl |M 3 | 2 - A, 2 \M,\\ (1.14) 

16tt 2 | M 2 s = X T - Qgl |M 2 | 2 - \ g\ \M X \\ (1.15) 

16vr 2 |Ml 3 = 2X r -^ 1 2 |M 1 | 2 , (1.16) 

16vr 2 - = 3X b + X T - 6g 2 2 |M 2 | 2 - - g 2 , \M,\ 2 , (1.17) 

16tt 2 ^ = 3X, - 6g 2 2 \M 2 \ 2 - ^g{ \M,\ 2 , (1.18) 

with 

X t = 2\h t \ 2 (m 2 H2 + M? 3 + M? 3 ) + 2K| 2 , (1.19) 

X b = 2\h b \ 2 (m 2 Hi + M| a + M? 3 ) + 2\A b \ 2 , (1.20) 

X r = 2\h t \ 2 (m 2 Hl + M? a + M| 3 ) + 2|A T | 2 . (1.21) 

X t , Xf,, and X r are always positive, so they decrease the scalar masses as one runs downwards 
to the low scale. Therefore, the soft breaking parameters of the third generation are in general 
smaller than those of the first and second generation [unless one starts with very different values 
at the high scale]. 

The RG equations for m 2 H and rn? H are still a special case. Compared to those for the soft 
squark masses they do not get the large contributions proportional to |M 3 | 2 . Moreover, X t 
and X b enter with larger coefficients. This can cause m/jj or m? H2 to become negative near 
the electroweak scale [pi . Since this leads to a breakdown of the electroweak symmetry solely 



by quantum corrections this is called radiative electroweak symmetry breaking ||15|| . Here we 
note that owing to the large top quark mass one can expect an especially large effect from X t , 
favoring m 2 H2 < rn 2 H (unless h t ~ hb). 



The gaugino masses do not only enter Eqs. (|1.11|) to (1.18) but also the RG equations for the 



A parameters. Non-zero gaugino masses at Mx are therefore sufficient to create all the other 
soft SUSY breaking terms. On the other hand, if the gaugino masses vanished at Mx they 
would only be generated through 2-loop and higher order effects and hence be very small. 

Another (technical) remark seems appropriate at this point: For computing radiative correc- 
tions in supersymmetry, it is important to choose regularization and renormalization schemes 
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that preserve supersymmetry (and gauge symmetry, of course). Dimensional regularization 
(DREG), for instance, explicitly violates SUSY because the continuation of the number of 
spacetime dimensions to D = 4 — e introduces a mismatch between the number of gauge boson 
and gaugino degrees of freedom. A solution is to perform the momentum integrals in D = 4 — e 
dimensions while taking the vector index \x of the gauge boson fields over all four dimensions 
IB|, |17| , 18| . This is known as dimensional reduction (DRED) and nicely respects both gauge 



symmetry and SUSY, at least up to 2-loop order. We will come back to this when discussing 
super symmetric QCD corrections in Chapter 3. 



1.3 Models of SUSY Breaking 

Above we have introduced explicit SUSY breaking terms because we are ignorant of the funda- 
mental machanism that breaks supersymmerty. If SUSY is broken spontaneously there exists 
a Goldstone fermion called the goldstino. In global supersymmetry the goldstino is massless. 
In local supersymmetry (supergravity) the goldstino is "eaten" by the gravitino Q73/2) which in 
this way acquires a mass 1713/2 [0- This is called the super-Higgs mechanism and is completely 
analogous to the ordinary Higgs mechanism in gauge theories. 

Present models of spontaneously-broken low-energy supersymmetry assume that SUSY is bro- 
ken in a "hidden" or "secluded" sector which is completely neutral with respect to the SM 
gauge group. The information of SUSY breaking is then mediated to the "visible" sector, 
which contains the MSSM, by some mechanism. There are no renormalizable tree-level in- 
teractions between the hidden and visible sectors. Two scenarios have been studied in detail: 
gravity-mediated and gauge-mediated SUSY breaking. 

In gravity-mediated SUSY breaking, SUSY breaking is transmitted to the MSSM via grav- 
itational interactions [20]. The breakdown of the symmetry occurs at O(10 10 ) GeV or higher, 



and the gravitino gets a mass of the order of the electroweak scale. The simplest realization 
of such a framework is the minimal supergravity model (mSUGRA) [[I], ^TJ. In this approach, 
one assumes a universal gaugino mass M1/2, a universal scalar mass Mq, and a universal cubic 
coupling A at Mx- In addition, one just needs to specify tan/3 and the sign of /zo- RGEs are 
then used to derive the MSSM parameters at the electroweak scale. Since mSUGRA involves 
only five parameters (in addition to the 18 SM parameters) it is highly predictive and thus 
used for most experimental searches. However, one should keep in mind that it also highly 
restrictive. 



Gauge-mediated SUSY breaking (GMSB) |22j models involve a "secluded" sector where 
SUSY is broken and a "messenger" sector consisting of particles with SU(3)xSU(2)xU(l) 
quantum numbers. The messengers directly couple to the particles of the secluded sector. 
This generates a SUSY breaking spectrum in the messenger sector. Finally, SUSY breaking is 
mediated to the MSSM via virtual exchange of the messengers. A basic feature of such models 
is that SUSY is broken at much lower scales than in the gravity-mediated case, typically at 
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0(1O 4 -1O 5 ) GeV. Moreover, the gravitino gets a mass in the eV to keV range, and is therefore 
the LSP. This can be crucial for SUSY signatures at collider experiments because the next- 
to-lightest SUSY particle (NLSP) will eventually decay into its SM partner plus a gravitino. 
A long-lived x°-NLSP that decays outside the detector leads to the usual SUSY signature of 
large missing energy plus leptons and/or jets. If, in contrast, the decay Xi 193/2 occurs 
inside the detector SUSY events would in addition contain photons. The NLSP may, however, 
also be a charged particle e.g., fjjr. This would lead either to a long-lived charged particle or 
to SUSY signatures characterized by r-leptons. 

Since gauge interactions are flavour-blind one has universal boundary conditions in GMSB as 
in mSUGRA. The low-energy spectrum is determined by the mass of the messengers. Minimal 
GMSB is thus even more restrictive than mSUGRA. In the most general case, however, both 
supergravity and gauge-mediated effects may contribute to the breaking of supersymmetry. 

A more detailed discussion of SUSY breaking is beyond the scope of this thesis. For a 
thorough introduction to supergravity, see e.g. 0; for the phenomenology of mSUGRA, see 



e.g. [23, M. A review of gauge-mediated SUSY breaking is given in 



1.4 Electroweak Symmetry Breaking 

The scalar potential for the Higgs fields — including all soft-breaking terms — is 

VHi gg s = (H 2 + m 2 Hl ) \H Y \ 2 + (M 2 + m 2 H2 ) \H 2 \ 2 + {BuHxH* + h.c.) 

+ l(gl + gl) {\H,\ 2 - \H 2 \ 2 ) 2 + \ gl \H+H° 2 * + H° 2 H^\ 2 . (1.22) 

Here the terms proportional to |/i| 2 come from F-terms, the quartic interactions come from 
.D-terms, and the terms proportional to m 2 H , m 2 H and B/i come from soft SUSY breaking, see 
Eq. 0. 

Fist, we use SU(2) L gauge transformations to rotate away any VEV of one of the charged Higgs 
fields e.g., (Hi) = 0. Then <9Vmggs/ 'dHi = requires (H 2 ) = as well. This is good because 
we are now sure that electric charge is conserved in the Higgs sector. We shall thus ignore the 
charged components in (1.22) when minimizing the potential []. 

Next, we choose Bfi, the only term in (1.22) that depends on complex phases, to be real and 
positive. This can be done through a redefinition of the phases of Hi and H 2 . Then (H®) 
and (H 2 ) are also real. This means that CP is not spontaneously broken by the Higgs scalar 
potential. The Higgs mass eigenstates are thus also eigenstates of CP. 

The scalar potential has a lokal minimum other than its origin if 

(i^r+< 1 )(H 2 +< 2 )<(^ ) 2 - (i-23) 

4 To be completely honest, this is a simplifying view: charge might still be broken in the absolute minimum 
of the full scalar potential due to some non-zero sfermion VEVs. 
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However, this is not enough; Vmggs must also be bounded from below. This is the case if 

(H 2 + m 2 Hl ) + (|/i| 2 + m 2 H2 )> 2Bfi. (1.24) 

These two conditions can only be satisfied simultaneously if 7^ m# 2 — implying that in 
the MSSM electroweak symmetry breaking is not possible without first breaking SUSY! 

We are finally ready to minimize the Higgs potential by solving dVuigg S /dH® = dVm ggfi /dH% = 0. 
The VEVs at the minimum of the potential, Vi = (Hf) and v 2 = {H®), are related by 

v 2 = v\ + v 2 = 2m 2 z /(g 2 1 + g 2 2 ) » (174 GeV) 2 , (1.25) 

so that only the ratio of the two remains a free parameter. Defining 

tan = v 2 /v 1 (1.26) 

the minimalization conditions are given by 

|/i| 2 + m Hl = —B\i tan — \ m 2 z cos 2/3, (1-27) 
|^| 2 + m 2 H2 = -Bficot(3 + ±m|cos2/3. (1.28) 

The two complex scalar Higgs doublets consist of eight degrees of freedom, three of which are 
eaten by the longitudinal modes of the Z and W bosons. The remaining five physical degrees of 
freedom form a neutral CP-odd, two neutral CP-even and two charged Higgs bosons denoted 
by A , h°, H°, and H ± , respectively. They are given by 

( G °\-,/?( Sin ^ 008 P\ (^ H ° 2 A n 001 

\A°) ~ VZ {-cos (3 sin (3 ) {^H?))' [ ^ 

(#+) = (-cos? Unp) (hI*) (L30) 
with G° and G + the would-be Nambu-Goldstone bosons, and 

cosa sina \ f^R(H 2 ) — v 2 \ , , 

\H°J ~ V ^-sina cosa J \&(H?) - vj 1 J 

with the Higgs mixing angle a. ^St(H) denotes the real and $s(H) the imaginary part of H.] 
The mass eigenvalues at tree level are 

m 2 A = 2Bfi/sm2[3, (1.32) 
m 2 H ± = m 2 A + mly, (1.33) 

m l°,H — \ ( m A + m l T \J ( m A + m 2 z ) 2 - Am 2 A m 2 z cos 2~Pj . (1-34) 
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We can therefore take tua and tan (3 as the free parameters of the MSSM Higgs sector. From 
(1.34) it follows that 

rriho < mm(m z ,mA) |cos2/3| (1.35) 

at tree level. However, the Higgs masses and mixing angle are subject to large radiative 
corrections. In this thesis we use the formulae of Ref. |26[ for the (h°, H°) system and Ref. 



for the mass of H (]. Assuming that the masses of the sparticles in the loops do not exceed 
1 TeV one obtains an upper bound for the lightest Higgs mass of 

m h o < 130 GeV. (1.36) 

On the other hand, A , H°, and H can be almost arbitrarily heavy. If m A 3> rriz they become 
nearly degenerate and decouple from the low-energy regime. In this case h° is very difficult 



to distinguish from a Standard Model Higgs boson. Figure |1.1| shows the MSSM Higgs boson 
masses as a function of uia for tan /3 — 3 and 30. For the radiative corrections we have used 
Mq = 400 GeV, M = 350 GeV, M b = 450 GeV, A t = A b = -300 GeV, and /i = 400 GeV. 
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Figure 1.1: MSSM Higgs boson masses as a function of for tan (3 = 3 and 30 (Mq = 400GeV, 
Mjj = 350 GeV, M D = 450 GeV, A t = A b = -300 GeV, and /i = 400 GeV). 

Up-type quarks get masses proportional to v<i = v sin/?; the masses of down-type quarks 
and electron-type leptons are proportional to v± = v cos (3. At tree level the quark and lepton 
masses are therefore related to the respective Yukawa couplings by 



K 



gm d 



m w sin (3 ' 



V2 m w cos (3 ' 



V2 m w cos (3 



;i.37) 



This is why we neglected the Yukawa couplings of the first and second generation, but not 
those of t, b, and r. Obviously, h t is significant due to the large top quark mass. Moreover, 

5 Notice that |2(| |27j have the opposite sign convention for the parameter fj,. 
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h t : hb : h T = m t : tan /3 : m T tan /3 from ( 1.37|) , so also the bottom and tau Yukawa couplings 
are important if tan (3 is large. In fact, certain models predict the unification of top and bottom 
(or top, bottom, and tau) Yukawa couplings for tan/5 ~ m t jm^. One can, moreover, require 
that none of the Yukawa couplings become nonpertubativly large. This gives rough bounds of 
1.2 < tan/? < 65. 



1.5 Sparticle Mixing and Mass Eigenstates 

Once SU(2)l x U(l)y symmetry is broken, fields with the same SU(3) C x U(l) em quantum 
numbers (and, of course, R-parity and spin) can mix with each other. In the Standard Model, 
B° and W l mix to 7, Z°, and W . Also the Dirac masses of quarks and leptons can be 
understood as such mixing terms. In the MSSM, this mixing also effects squarks, sleptons, 
Higgs bosons, as well as gauginos and higgsinos. The lone exception is the gluino, being the 
only colour octet fermion in the model. 

Indeed, masses and mixings of sparticles are of crucial importance both theoretically and ex- 
perimentally: i) they determine the properties of the sparticles searched for and ii) they are 
directly related to the question of how SUSY is broken. (Notice that the main questionmarks 
in the MSSM come from £ so f t while the couplings and all but one mass term in £susy are linked 
to well known parameters of the SM!) 



1.5.1 Squarks 

In the most general case, the squark mass eigenstates are obtained by diagonalizing two 6x6 
squark mass-squared matrices — one for up-type and one for down-type squarks. However, 
mixing between squarks of different generations can cause severe problems due to too large loop 
contributions to flavour changing neutral current (FCNC) processes [^HJ . 
Ignoring intergenerational mixing for the moment, the two general 6x6 squark mass-squared 
matrices decompose into a series of 2 x 2 matrices, each of which describes squarks of a specific 
flavour [p9|: 



M- 



with 



m~ L a q m q 
a q m q m\ R 



(JZ 



'ht I "'"i: 



771 



;i.38) 



777 



m 



<iL 
2 

a„ 



Ml + m| cos 2/3 (i| L — e q sin 2 9 W ) + m 



Mffj + e q m| cos 2(3 sin 2 W + m 
A q — fi {cot f3, tan/?} 



;i.39) 

;i.4o) 
;i.4i) 
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for {up, down} type squarks, respectively. e q and If are the electric charge and the third 
component of the weak isospin of the squark q, and m q is the mass of the partner quark. Mq, 
Mjy, and are soft SUSY breaking masses, and A q are trilinear couplings as in ( |1.8|) . Family 
indices have been neglected. 

The off-diagonal elements of M 2 ^ are proportional to the mass of the corresponding quark. In 
the case of the 1st and 2nd generation q L and q R are therefore to a good approximation also 
the mass eigenstates. Indeed, experimental results on K° — K° and D° — D° mixing suggest 
that the masses of up, down, charm, and strange squarks are highly degenerated and mixing 
can be neglected [PSfl . We can therefore assume 



mu L = = = mg R = m- CL = m- CR = m- SL = m- SR . (1.42) 

This also justifies the above assumption that intergenerational mixing can be neglected. 

However, this does not hold for the third generation: Stops are expected to be highly mixed 
due to the large top quark mass, and for sbottoms mixing effects can be important if tan (3 is 
large. In the following we discuss this left-right mixing for q = t,b: 



According to Eq. ( |1.38| ) M. 2 q is diagonalized by a unitary matrix 1Z^. Assuming that CP violating 



phases only occur in the CKM matrix, we choose lZ q to be real. The weak eigenstates q L and 
q R are thus related to their mass eigenstates q x and q 2 by 

(?)=WH *>'=( C ° S ^ (1-43) 
with 9q the squark mixing angle. The mass eigenvalues are given by 



91,2 



\ (™& + ™\ R T VK ~ m V 2 + Aa l m t) ■ ( L44 ) 

By convention, we choose q ± to be the lighter mass eigenstate. Notice, that m qi < m qL R < m q , 2 . 
For the mixing angle 9 q we require < 9 q < it. We thus have 



Moreover, | cos6>g| > if m qL < m qR , and | cos 6^1 < 4g if m qR < m qL . 



Inverting Eqs. (|1.38| ) - (|1.44j ) one can calculate the underlying soft-breaking parameters from 
the physical squark masses and mixing angles: 



• cos 2 6i + m| 2 sin 2 9i — m| cos 2(3 (| — | sin 2 9w) 



M% = m 

M^j = m\ i sin 2 6$ + mf cos 2 — § m l cos 2/? sin 2 9w — rn 



rn 



1 1 



in the stop sector and 

Mq = m 2 x cos 2 6- b + m~ b2 sin 2 9 b + m| cos 2(3 (| - § sin 2 9 



Ml 



m 



• sin 2 9 b + m 2 2 cos 2 9 b + | m| cos 2/5 sin 2 9 



w 



w , 



m 



b- 



1.46) 

;i.47) 
;i.48) 

1.49) 
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in the sbottom sector. Finally, the off-diagonal element of the squark mass matrix is 

a q m q = ~ (m? - m~ 2 ) sin 20, . (1.50) 

Notice that the parameter Mq enters both the stop and the sbottom mass matrices. Therefore, 
(1.46) and (1-48) imply the following sum rule at tree level: 

m 2 ^ cos 2(3 = m| cos 2 % + mf 2 sin 2 9i — mj- cos 2 6 b — mf^ sin 2 6 b — m 2 + m 2 . (1.51) 

This shows that if tan (3 and five of the six quantities , m~ t2 , m bi , m~ b2 , % 9- b are known, the 
sixth can be predicted. 

In Sect. |1.2| we have learnt that the soft SUSY breaking parameters (and hence rrtq L and 
rriq R ) of the 3rd generation are most likely different from those of the other generations. Left- 
right mixing induces an additional mass splitting, rendering q x lighter and q 2 heavier than q L R . 
Notice that both effects are induced by large Yukawa couplings. The mixing angles (and Yukawa 
couplings) also enter the stop and sbottom couplings to other particles. Moreover, stops — and 
for large tan (3 also sbottoms — give important contributions to radiative corrections in the 
Higgs sector (think e.g. of radiative symmetry breaking!). 

There are some experimental constraints on the mass splitting in the stop and sbottom sec- 
tors from BR(6 — > 57) and Sp. However, these constraints are much weaker than the above 
mentioned constraints from K° — K° and D° — D° mixing. 

All together, we now understand that the phenomenology of stops and sbottoms can be very 
different from that of the other squarks. 



1.5.2 Sleptons 

The mass matrix of the charged sleptons is completely analogous to that of the squarks: 

K = L i L I2 , l = S, fi, r, (1.52) 



with 



M 2 L -m| cos 2/3 (± -sin 2 9 W ) + m 2 , (1.53) 



mj R = Ml -m 2 z cos2(3 sin 2 6 W + m 2 e , (1.54) 
ai = Ai — pta.n/3. (1.55) 

From RG evolution one expects M| < M| and hence m^ fl < rn^. For selectrons and smuons, 
the "left" and "right" states (eL,R and Pl,r) are also the mass eigenstates. For staus, however, 
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analogous arguments apply as for the sbottoms. If tan (3 is large enough, r L and t r therefore 
mix to mass eigenstates 

f x = cos^Tl + sin 6^7^, (1.56) 
f 2 = cos^ff^ — smOff L (1-57) 

with Of the mixing angle. All formulae given for squark mixing are also applicable in this case. 

Let us now turn to the sneutrinos. In the case of massless neutrinos, there is only one 
sneutrino, u L , with a mass 

mo L = Ml + \ m| cos 2(3 (1.58) 

for each generation. If, however, neutrinos have a mass (as suggested by the Superkamiokande 
experiment |JIJ) there exists also vr. We assume that the neutrino masses are generated by 
the see-saw mechanism 

(•*■ *) ( z z: ) ft) 

where iV is a heavy singlet field ||. tum is the VjV^ Majorana mass and mo the vlN Dirac 
mass. Mm, the Majorana mass of N, is very large. In fact, it can be as large as the GUT scale. 
Therefore, we may expect that vr is also very heavy and does not contribute to low-energy 
phenomenology. 



1.5.3 Charginos 

The supersymmetric partners of the W ± and the H ± mix to mass eigenstates called charginos 
xf (i = 1, 2) which are four-component Dirac fermions. In order to deduce the properties of 
the latter we start with the basis 

^ = (-zA + , Vk) , r = (-»A- iJ> 2 Hl ) , (1.60) 

where A ± = 4g (A 1 =F ^A 2 ), tp}j 2 , and ip 2 Hl denote the two-component spinor fields of W ± ) H + , 
and H~, respectively. The mass terms of the lagrangian of the charged gaugino-higgsino system 



can then be written as 0, 31 



withQ 

v — / 

m w cos (3 n 
3 Here and in the following M = M 2 



x= ( M ^2m w sm(3\ (162) 
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The mass matrix X, Eq. (|1.62| ), is diagonalized by two 2x2 unitary matrices U and V: 

WXV- 1 = M c (1.63) 

with Ai c the diagonal mass matrix. With the rotations 

xt = V tj ^, X 7 = U ij ^, 1,3 = 1,2, (1.64) 

the mass eigenstates in Dirac notation are given by 

Xt=(€), (1-65) 



We take xl to be the lighter chargino per definition. Moreover, assuming CP conservation, 



we choose a phase convention in which U and V are real. Following Ref. |3T[ the mass eigenvalues 
then are 



M n = UXV 







;i.66) 



with r/ i = ±1 the sign of the eigenvalue, and m-+ = |(M c )jj| the chargino mass. According to 
Eq. ( jTBg ) we have (i = 1,2) 



(M ( 



Oii — 2 



(M -fi) 2 + 2 m 2 w (1 + sin 2/3) t J (M + fi) 2 + 2m 2 w (I - sin 2(3) 



1.67) 



The larger mass eigenvalue is always positive, rj 2 = 1. The smaller eigenvalue is positive if 
M-fi — m 2 ^ sin2/? < 0, and negative otherwise. 

The elements Uij and Vij of the diagonalizing matrices can also be directly expressed by the 
SUSY parameters M, /i, and tan/?: 



0i L , M 2 -/i 2 -2mH,cos2/3 

f/l2 = f/21 = 71V 1 + w 



6 2 I M 2 - fi 2 -2m w cos2(3 

u ™ = - Ull = 7 = 2V w 



v 21 = -V 12 = -^\ll 



6 3 I M 2 — jj, 2 + 2 m w cos 2(3 



s/2 



W 



9 4 I M 2 — jj, 2 + 2 m w cos 2(3 

V22 = Vll = T2T w 



;i.68) 
;i.69) 
;i.ro) 
;i.7i) 



with 



W = J (M 2 + /i 2 + 2m 2 v ) 2 — 4 (M-fi — sin2/3) 2 



;i.72) 
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and the sign factors 0j, i — 1 . . . 4 

{9i, 6*2, # 3 , 64} = 



{1, e fl , !} • • • tan/? > 1 
{£ fl , 1, 1, £ A } ... tan/3 < 1 



where 



e A = sign(M sin (3 + n cos /?), e s = sign(M cos (3 + /i sin /3). 



(1.73) 



1.74) 



In what follows we will always use the convention (|1.63|) ||31|| . However, one can also choose 
U and V such that M. c has only non-negative entries 0. This is for example achieved by 
solving the eigenvalue problem for X T X, 



Mr = diaff(m 2 - + , m 2 ~+) = VX L X V 



TvtM 



with 



V 



xi 1 xi 



cos</) 1 sinipj 



;i.75) 



;i.76) 



and the matrix U given by 



U — ~~ri — V X 



T 



M 



The mass eigenvalues then are 



;i.77) 



m lt = 5 



M 2 + /i 2 + 2m^ 



=F J(M 2 -^) 2 + Am 2 w cos 2 2f3 + Am 2 w (M 2 + /x 2 + 2 M-/i sin 2/3) 



1.78) 



This can be a practical approach when treating the problem numerically (the analytic ex- 
pressions for U and V become rather complicated) and is especially useful in case of complex 
parameters, see e.g. [32]. However, one has to take care that <p 1 and <fi 2 are always in the correct 
sector of the unitcircle e.g., G [0,7r). 



1.5.4 Neutralinos 

The neutral gauginos and higgsinos also mix. Their mass eigenstates are the neutralinos Xi\ 
z = 1 ... 4. In general, both weak and mass eigenstates are Majorana fermions; however, if two 
neutralinos are degenerated in mass they can combine to a Dirac spinor. 
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In this thesis we choose the basis 



with 



= ^Hj sin P - V'ia COS P, = 4>Hi cos P - ^H 2 Sin P, 



1.79) 



;i.8o) 



and A 7 , A z , ipjj , 1^% the two-component spinors of the photino 7, zino Z°, and the neutral 
higgsinos Hf, H®, respectively. The mass terms of the neutral gaugino-higgsino system can 
then be written as 0, ^3] 



r 



■!(v>°) T y# ) + h.c. 



:i.si^ 



with 



y 



( M'c 2 w + Ms 2 w (M-M')s w c w 
(M-M')s w c w M's 2 w + Mc 2 w 
m 7 



rrir 



^ s 2f3 
~L LC 2(3 



\ 




;i.82) 



\ 

where we have used the abbreviations s 
M and M' are the SU(2) and U(l) gaugino masses, M = M 2 , M' = M 1 (we will stick to this 
notation in the following). 



sin 6^, c w = cos Ow, $2/3 = sin 2/5, c 2/ 3 = cos 2f3. 



The mass matrix Y is diagonalized by a 4 x 4 unitary matrix N, 
N*Y N- 1 = M, 



:i.83i 



with M. N the diagonal mass matrix. The mass eigenstates in two-component notation then 
are 





13 rj ' 



1,3 



and in four-component notation 



\Xi J 



1...4. 



;i.84) 



;i.85) 



As in the charged gaugino-higgsino sector one can choose N such that Ai N has no negative 
entries by solving the squared eigenvalue problem, 



diag(m| ? , ... , mjo) = NY^Y N~ 



1.86) 



However, assuming CP conservation, we find it convenient to allow for negative mass eigenval- 
ues. We can then choose a phase convention in which iV is real and orthogonal: 



M 



N 



NY N = diag(e 1 m^o, e 2 m^o } £ 3 m^o, ^m^o) 



"XV 



;i.87) 
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with £j = ±1 the sign of the eigenvalue and m^o = \(M.n)u\ the mass of the neutralino xl 
(0 < m^o < rn^o < m^(> < m^o by convention). 

In the case M <C the two lightest neutralinos are dominated by their gaugino compo- 
nents. In particular, Xi ~ B and Xi ~ W^ 3 - Similarly, Xi is mostly a charged WMno. For the 
masses one finds very roughly 2m^o ~ m^o ~ m-+ ~ M. In the opposite case <C M), x° 2 
and are mostly higgsinos with masses close to Finally, for M ~ the gaugino and 
higgsino states are strongly mixed. 



1.6 Interaction Lagrangian 

In this section we list the relevant parts of the Lagrangian and the Feynman rules for squark 
(and related) interactions needed in this thesis. All interactions are given in the (q L , q R ) as well 
as in the (qi,q 2 ) basis. We concentrate on the 3rd generation (q = t,b; q — t,b). However, all 
formule also apply for the 1st and 2nd generation, provided one inserts the (super) CKM-Matrix 
elements in a proper way. For the 3rd generation we take (CKM) 33 = (sCKM) 33 = 1. The 
corresponding expressions for (s)lepton interactions can be derived by the obvious replacements 
q — > £, t — > v, b — > £~, q — > £, e q — > e e , etc. 



1.6.1 Quark — Quark — Gauge Boson 

£ qqi = -eegA^q^q (1.88) 
£ qq z = ^— Z^qY {(i| L - e q sin 2 6 W ) P L - e q sin 2 6 W P R } q 

COS U\y 

= --77T Z » ( C iL P l + C qR P R ) q (1.89) 

COS U\y 

C qL ,R ■= IIl,R~ e 9 sin2 °W (1-90) 



■-qqW 



"71 * ^ Pl 6 + w * 1 r Pl t] (L91) 



C qq9 = -g s T r a s G;q r Yq s (1-92) 
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1.6.2 Squark — Squark — Gauge Boson 

(a) squark - squark - photon 

Cqrr = iee q Ap(? L irq L + fa& t q R ) 

= iee q A^n\{R)^n\ 2 n%)~q)^~ qi 

= iee q 8 ij A ll ^& l q i (1.93) 



(b) squark - squark - Z° 

C ^ z = T^r Z » ^ ft & h + c <* & & ^ = T^r °a Z » % & * (L94) 

COS uw cos "w 
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c ij '■— C qL 7l il Tlj l + C qR IZ i2 'R.j2 



I\ L cos 2 9g — e g sin 2 6 W —\ J| L sin 26 q \ 
- \ J L sin 26 q II L sin 2 6 q - e q sin 2 6 W ) .. 



(1.95) 



(c) squark - squark - W ± 

Aw = ^(w^t* L dn L + w-b* R dH R ) 

= ^(^5 1 W7??6 i + ^J 1 W'-6;?f i ) (1.96) 

(d) squark - squark - gluon 

The corresponding Feynman rules we obtain from 

A&B = A (d^B) - (d^A) B -> q* q t = % (k t + k^ (1.98) 
where k i and kj are the four-momenta of q t and q~j in direction of the charge flow. 



k / * 
x 



7 vAAAAAAAA*.^ — iee g (j? + fc) M 5jj 



Z vAAAAAAAA*.' c (p _)_ MM 

cos 0w ' 
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k , 



w ± ^wwv< -j. k\ x k% (p + ky 



/ 

g ^mmrncex^ -i g s T r a s (p + k)^ 8, 



1.6.3 Quark — Quark — Higgs Boson 



C qqH = slh qq + s q 2 H qq + slA q^q 



+ H + i(siPr+S b A P l 



R) b + H-b(s\P L + s\P R )t 



1.99) 



with 



s = - 



gmt 



b 2 



rt COS CX « 

2 mw snip ' 

- 9 gmt . s sin a, 

2 raw sin /3 ' 



s| = i^-cot/3, 

•5 2 ^ ' 



i>4 



gm t 



mw 



V2 



cot (5 



b 2 



9 g sin a, 

2 m)4/ cos p ' 

7i ^ 77 COS Gi ^ 

2 mw cos p 



„6 
A 4 



rati/ 



tan/3. 



1.100) 
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1.6.4 Squark — Squark — Higgs Boson 

Defining H k = {h°, H°, A , H ± } we can write the Higgs interaction with squarks in the general 
form 

Cm = H k (t, 4*) Gk (|) = (G fe % H k (1-101) 

where a and j3 are flavor indices. For k — 1, 2, 3 we have of course a = [3; in case of k = 4 we 
have G 4 = Gf = (Gf) T . G fc and G k are related by 

Gl = K*Gl{1Vy, (A; = 1,2,3) (1.102) 
G* = K l G 4 (7^) T , G\ = U b (G 4 ) T (TZY = (G|) T (1.103) 



~0 
^3 



n a g? (npy- 
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The individual contributions are: 



(a) squark - squark - h° 



£qqh° = 9 { [ 



^(i-§sin 2 Msin(a+/3) 



raw sin 



COS OL 



h°t* L t L 



+ 



2 ra z 



sin 2 9 W sm(a+/3) 



cos a 



raw sin 

,or* r . 7 or* 



3 cos 0i4^ 

(A cos a + ix sin a) (h%t L + fe ^) 



^RpR 



^(|-|sin 2 Msin(a+/5)- 



+ 



1 _mg_ 

3 cos 0w 



ra t 
2mw cos /3 



sin 2 9w sin(a+/3) 



raw sin 



raw cos /3 



sin a 



sin a 



(A b sin a + cos a) (h°b* R b L + h°b* L b R ) j 



and thus 



"g '"q "I _ 



-71 ^o* 

- a K { } + /* { 4: } ) ^ ^ ^ - ^ { } 



(1.104) 



;i.io5) 



for | } type squarks respectively. We use the abbreviations c w = cos9 w , s q = sin a, 

c a = cos a, s a+/ 3 = sin(aH-/3), = I$L—e q sin 2 0^, and = e q sin 2 6*^. a is the mixing angle 

in the CP even neutral Higgs boson sector. h q are the Yukawa couplings: 



ht 



gm t 



72 mw sin (3 ' 



gm b 



\^2m w cosf3 



1.106) 



(b) squark - squark - H° 

C m o = -g {[-^-(I-|sin 2 Mcos(a+/5) + 
+ 



mw sin/3 



sin a 



H°t* L t L 



■ sin a 



+ 



raw sin 

— (A sin a - fi cos a) (#°^~ L + H°i* L i R ) 



H °t*RpR 



2raw sin 



-^f- (± - ± sin 2 W ) cos(a+/3) ^ 

cos 0jy V 2 3 " / V " / mvi/ cos /3 



lc^sin 2 ^cos(a+/5) 



• cos a 



cos a 



H%b R 



H°~b* L b L 



+ 



rab 



Iraw cos 



raw sin 

(A b cos a - (j, sin a) (H°b* R b L + H%* L b R ) } 



(1.107) 
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gm. 



C qL c a+ , - V2 m q h q { } -^M^iZ}-^}) 



- H 



9 m Z C C n 



\[2m q h q 



1.108) 



Notice that G\ can be obtained from G\ by the replacement a — > a + 7r/2, i.e. c a — > s a and 



(c) squark - squark - 1° 



-■qqA 



2m 



w 



m t (A t cot (3 + fi) (A%t L - A%t R ) 

+ m h (l 6 tan/3 + //) (A%b L - A°^6 fl ) 



(1.109) 



° 2 mw 



-A 



cot/3 
\^ q \ tan/3 



cot/3 
1 tan/3 





(1.110) 



In this particular case we have G| = G§ because Gf ist an off-diagonal and 1Z? is an unitary- 
matrix. 



(d) squark - squark - H 



± 



qqH± 



V2 



m 



w 



(m 2 b tan (3 + m\ cot (3 - m 2 ^ sin 2/5) + t* b L + ^™ 



+ m t (A t cot /3 + fi) H + t* L b R + m 6 (A b tan /3 + //) 6 L 1 + h.c. (1.111) 



ml tan (3 + m 2 cot (3 — m^ sin 2/5 m& (A\ tan (3 + /i) 
G A = ^J^\ ' ] (1.112) 

^ m t (A t cot (3 + 
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1.6.5 Quark — Squark — Chargino 



gt{-U l3 P R + Y t V 2j P L ) xf b L + gi(Y b U 2j P R ) j$ b R 
+ g b (- Vy P R + Y b U 2j P L ) xf h + gb (Y t V 2j P R ) xf h 
+ g l] (-C/ y P L + Y t V 2] P R ) tbl + g l] (Y b U 2j P L ) t b R 



~,+ 

■3 

3 



+ g xT (~Vi3 Pl + Yb U 2j P R ) bt* L + g l] c (Y t V 2j P L ) b t* R 

gHit j P R + klP L )tfb i + gb(4 j P R + kj j P L )xl c t i 
+ 9xU i l P L + klP R )tbt + g^ c (£l j P L + ki j P R )bt* 



(1.113) 



The (Ji-q'-Xj couplings £fj and kL can be written as 



(1.114) 



with 

1 {Y,V j2 )' U > UtW 

and the Yukawa factors Y q = h q /g, i.e. 

m t m b 



Y, 



m w sin (3 



Y b 



m w cos (3 



(1.115) 



(1.116) 



U and V are the 2x2 unitary matrices diagonalizing the charged gaugino-higgsino mass matrix 



as defined in Sect. 1.5.3 




i9 (4 Pr + KjPlj 




ig (4 P R + k% P L 
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i9 (4 Pl + ^Pr) 



ig (4 P L + fc| P fl ) 



1.6.6 Quark — Squark — Neutralino 

£mx° = 9q(fl k PR + h Lk p L)xlh + gq(h% k p R + f q Rk p L )x° k q R + h. c . 

= g^{o q ik P R + b q ik P L )xlq i + g xl (a 9 i k p l + b lk p r) <1 V* (1.117) 
The (ji-q-Xk couplings a? fc and are given by 



E *t ^L, 4 = E *t ^1 (1.118) 



with 



and 



/Ife = -V2e g sin%AA fel - V2 (I q 3L - e q sm 2 9 w )£^, (1.120) 

/l fc = -V2e q sm6 w {tSLn6 w N k2 -N kl ), (1.121) 

^ = K t (iV fc3 sin/3 - iV fc4 cos/?) = h Rk , (1.122) 

/i 6 Lfc = -y 6 (7V fc3 cos^ + Ar fc4 sin^) = h b Rk . (1.123) 



iV is the 4x4 unitary matrix diagonalizing the neutral gaugino-higgsino mass matrix, see 
Sect. Ol. 
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1.6.7 Interactions with Gluinos 

(a) quark - squark - gluino 



£ - 



-V2g s T r a s [(q r P R ~g a q LjS - q r P R g a q R J + Q a P L q r q* L , s ~ t P L q R ? R , S ) 
-V2g s T r a s \q r (niP R - TZf 2 P L ) g a ~ %s + ~g a (72«P L - K%P R ) q r 



(1.124) 



Note the relative minus sign between the terms with q L and q R : This is due to the facts that 
q R are colour anti-triplets and the anti-colour generator is — T at . 
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(b) gluon - gluino - gluino 



r .. -lil f . n a a h V o c 



(1.125) 



Owing to the Majorana nature of the gluino one must multiply by 2 to obtain the Feynman 
rule (or add the graph with g <-> g)\ 



/ 1i 

s 



-V2ig s T«(K! 1 P L -'Rl 2 P R ) 



Qi \ 

s 



-V2ig s T r %{niP R -T^ 2 P L ) 



g r mmmm 




9 s fabc 1^ 



1.6.8 Squark — Squark — Gauge Boson — Gauge Boson 



-<j<j77 — e 2 e q 1 M 1 M (q* L q L + q* R q R ) — e 2 e 2 q 5^ 1 M 1 M q* q i 



1.126) 



9' 



qqZZ 



COS 2 9 W 



9 



cos 2 9 W 

9 2 
cos 2 9w 



Z,Z»(C q 2 L q* L q L + C q 2 R q R q R ) 

z, z" (c q 2 L nin% + c 2 qR n%n%) q* $ 



z i:j Z^ Z» q* q i 



(1.127) 
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w 



* Qi 



Z •n" 



w 



/ 9j 



/3 



2 



7">9 -T-.9 



»3 7? 9 7? 9 a 
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2 iee q g s T a 8^ g^ 



V2tgg s T a Tli'Jlf 1 g^ 



1.6.9 Four— Squark Interaction 



The four-squark interaction stems from the D-terms in the scalar potential: 



y 2 



l\2 



D a D a + D i D i + 



+ 



(1.138) 



with a = 1, ... 8 and i — 1,2, 3. For our calculations we just need the SU(3) part \ D a D a which 



leads to: 

£qqqq = ~\ dl ^rs^tu (?L*r Q.L,s ~ Qli,r Q.R,s) {o.L,t Q.L,u ~ Qjt,t Q.R,v) 



(1.139) 



where a and (3 are flavor indices. Again a relative minus sign occurs between the q\q^ and 
q* R q R terms because the anti-colour generator is —T a * (see interactions with gluinos). In (q^ q 2 ) 
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notation (|1.139|) reads 



qqqq 



1 2 rpa rpa (fpc'jya 

2 "s 1 rs 1 tu V'^il'^jl 



1j,r %s 



(K 



kl'Hl 



n P npfl \ -13* ~3 



_1 2 rpa rpa qa qfl ~a* ~a ~/3* ~/3 
2 iJs J rs J tu °ij °kl Hj.r 1i,s %t Hk.v 



1.140) 



with 



S 



cos 28g — sin 29g 
— sin 20g — cos 29 g 



1.141) 



~/3 v 



9j 



rpa rpa c« CP I rpa rpa ca ca r 
1 rs 1 tu a ij a kl ~r 1 ru 1 ts a il °fcj°< 



a/3 



no sum over a. 



Chapter 2 

Squark Production at Lepton Colliders 



The next generation of high energy colliders (LHC, upgraded Tevatron, e + e~ linear colliders, 
colliders) will explore the TeV mass range testing the concept of low energy supersym- 
metry. While hadron colliders are well designed for broad searches, it is commonly expected 
that for a precise determination of the underlying SUSY parameters a lepton collider will be 
necessary. 

In this capter we discuss the pair production of squarks in e + e~ and fi + jjT collisions. For 
squark production at hadron colliders, see e.g. [^3| ^4], |34}] . 



2.1 Cross Sections for e + e Colliders 



The process e + e — > q { qj proceeds via 7 and Z exchange, see Fig. [2.1| a. For unpolarized 
beams the cross section at tree level is given by [[35], [J7[ 



tree V) 

U - A 



e q 5ij + 



1K r 4 4 



zz 



e q V e C ij S Jj n 



-w 



cos 2 Qw, and 



(2.1) 



where y/s is the center-of-mass energy, = [(s — m|. — m 2 ^) 2 — 4m? m| p 2 , s 2 ^ 



sin 2 fl 



w, 



D 



zz 



m 



2 \2 

z) 



7171 



2 ' 



D 



s (s 



rrtr 



7 z 



[s — m|) 2 + V z m z 



(2.2) 



v e = 4sin 2 ^vK — 1 and a e = — 1 are the vector and axial- vector couplings of electrons to Z 
bosons, and is the q^Z coupling as defined in (|1.95|) . The first term of Eq. (|2.1|) comes 
from pure 7 exchange, the second from pure Z exchange, and the third one is due to the 7— J? 
interference. Notice that the associated production of q x and q 2 [i 7^ j in ( |2.1|) 1 proceeds only via 
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Z exchange. The angular distribution shows the typical sin 2 $ shape, where $ is the scattering 
angle of q { : 



da 



tree 



d cos $ 4 



sin d a 



tree 



(2.3) 



The j—Z interference leads to a characteristic dependence of the cross section on the squark 
mixing angle In case of q 1 q 1 production the cross section is maximal for cos 6g = 1, 

i.e. q x = q L , and minimal for 



cos Of 



sin 2 6 



w 



l 3L 



1+1- 



4(t> 2 + a 



cos 2 6 W 



(2.4) 



In case of q 2 q 2 production the cross section is maximal for cos 6g = and minimal for 1 — 
cos 2 8q\ m i n . The cross section of q x q 2 production Q is maximal for maximal squark mixing 
(cos9q = ±l/v2) and vanishes in case of no mixing. Notice, however, that the sign of cosOq 
cannot be determined from cross-section measurements because the latter depends only on 
cos 2 9q. 

In case of a polarized e~ beam the total cross section reads: 



a(V-) 



7ra 2 K^ 



,2,$.. 



Q„2 2 
4 °w^W 



aV)D 



+ 



16 s^/Cyj/ 



[vl + a 2 ) 



2v e a e V- 



D 



zz 



(2.5) 



where V- is the degree of polarization, V- G [—1, 1]. V- = —0.8, for instance, means that 80% 
of the electrons are left-polarized while the remaining 20% are unpolarized. Thus, a(— 0.8) = 
0.9 a L + 0.1 a R with a L and cr R the cross sections for pure left-polarized (V- = —1) and pure 
right-polarized (V- = 1) e~ beams, respectively^. In general terms: 



[l-V-)<r L + {\ + V-)o 



R 



(2.6) 



If both the e+ and e beams are polarized we have 



e>8 ii {l-VJP+)- 



;i - V-V + ) - a e (V- - V + ] 



+ 



I J 



16 s^/C^/ 



{v 2 e + a 2 ) (1 - V-V + ) -2v e a e (7>_ - V + ) 



D 



D 



zz 



(2.7) 



Here P_ denotes the polarization factor of the e beam and V + that of the e + beam; V± 
{ — 1, 0, +1} for {left-polarized, unpolarized, right-polarized} e ± beams. 

1 Here and in the following q x q 2 means q^q 2 and e -9-i ^CSife) = a (QiQ2) + a (Q2<ii)- 
2 In the notation of fH this corresponds to £ = —0.9. 
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The radiative corrections we consider for the process e + e — ► g^g- are SUSY-QCD 
corrections to 0(a s ) and initial state radiation (ISR). Both have turned out to be significant. 

As for the SUSY-QCD corrections, they can be split into the standard QCD, i.e. gluonic, 
correction 6a 9 , the correction due to gluino exchange 5a 9 , and the correction due to squark 
exchange 5a q : 



a 



a 



r tree 



5a 9 + Sa 9 + 5a q . 



{21 



The corresponding Feynman diagrams are shown in Fig. |2.1| b-j. 
The gluon contribution Sa 9 has first been calculated in HO, 41]. It can be written as 
"4 a, 



5a 9 



3 7T 



-A, 



a 



tree 



(2.9) 



For the explicit form of Ay, see |p7[ . The SUSY contributions, 5a 9 and 8a q , were calculated 
in pB| , |3"7| ]. The inclusion of the gluino and squark exchanges in the on-shell renormalization 
scheme requires, however, a proper renormalization of the squark mixing angle. This problem 
has been solved in Ref. [37 which is thus the first complete treatement of supersymmetric 
QCD corrections to this process. We will discuss the subtleties of SUSY-QCD corrections in 



more detail in Chapter 3. Here we just note that within the scheme of [37] the counterterm 
for the squark mixing angle 86q is chosen such that it cancels the off-diagonal contribution of 
the squark wave-function corrections (i ^ k in Fig. |2.1| h, j). The squark contribution to the 
correction vanishes, 5a q = 0. The remaining correction due to gluino exchange is of the order 
of ±1% for an e + e~ collider with y/s < 200 GeV (LEP2). However, 5a 9 becomes important at 
higher energies as they are proposed for a Linear Collider: For large ^/s it can be up to —50% 
of the gluon correction. Moreover, 5a 9 does not scale with a tree and thus gives an additional 



dependence on the squark mixing angle. See [[37], [42] for more details and a numerical analysis 



The initial electron and positron may loose energy through photon emission. This is known 
as initial state radiation (ISR). The actual e + e~ annihilation thus takes place at the reduced 
center-of-mass energy s = s(l — x) with x the energy fraction carried away by the photon(s) 
Pi: 



a 



dx 



o 



+ 2/?) _ _ |) 



+ 7T \ 6 A, 



a (s(l - x)) 



where ar, is the cross section without ISR correction and (3 = 2a / '% (log s / ml — 1 
numerical evaluation of the ISR correction use the Monte Carlo routine PHOISR 



(2.10) 



For the 



Beamstrahlung also reduces the center-of-mass energy of the e + e~ collision. However, while 
ISR photons can in principle be dectected, beamstrahlung leads to an uncertainity in y/s. We 
thus leave this effect to Monte Carlo studies. 
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Figure 2.1: Feynman diagrams relevant for the process e + e~ — > q^y. (a) tree level, (b-f) gluon 
corrections, (g, h) gluino corrections, and (i, j) squark corrections. Lines that are not lettered 
in (b-j) are the same as in (a). 
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2.1.1 LEP2 



In 1995, after seven years of high precision measurements at the Z pole, the LEP accelerator 
was upgraded to its second stage LEP2 j45]. One of the main motivations was the search for a 
light Higgs boson and for SUSY particles. In particular, LEP2 covers an energy range that is 
hard to explore at LHC 



The LEP2 runs of 1995 to 1998 are listed in Table 371 [4TS]. The limits for SUSY searches 



presented at the ICHEP'98 conference in Vancouver are [47, 48]: 



m h o A o > 80 GeV, m H o > 90 GeV, and m H+ > 59 GeV, 



m,o > 28 GeV, 

m-+ > 90 (57) GeV for a gaugino-like xt an d tan/? = \/2 (40), 



m~+ > 62 GeV for a higgsino-like xt 



m e - B > 85 GeV, m~ u „ > 78 GeV, m fn > 72 GeV, 



m 



fi > 87 (84) GeV for B~ t = 0° (56°), and m~ bi > 87 (75) GeV for 6~ h = 0° (68°). 



With the data collected and analyzed since then, these limits have been improved by O(10 GeV) 



Date 




/ Cdt /exp. 


Nov. 95 & Oct. 97 


130 GeV 


6pb _1 




136 GeV 


epb^ 1 


July-Aug. 96 


161 GeV 


10 pb" 1 


Oct. -Nov. 96 


172 GeV 


10 pb" 1 


Aug.-Nov. 97 


183 GeV 


55 pb" 1 


May-Oct. 98 


189 GeV 


160 pb" 1 



Table 2.1: LEP2 runs in 1995 - 1998 



In 1999 and 2000 LEP2 will be operated at t/s ~ 200 GeV |gg|. Figure gja shows the 
iiti and bibi production cross sections at this energy for = 95 GeV as a function of 

cos9q (q = t,b). According to (|2.4| ) <r(titi) has its minimum at | cos %| ~ 0.5 and cr(6 1 & 1 ) at 
|cos^| ~ 0.35. In Fig. p72|b the corresponding ranges of t x and b\ production cross sections 



are shown as a function of the squark masses. For the SUSY-QCD corrections we have taken 



t-2 



m 



b 2 



m~ n 



500 GeV in both plots. 
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Figure 2.2: Cross sections for e + e — ► titi (full lines) and e + e — > 6161 (dashed lines) at 
x/s = 200 GeV ( function of cos 9q for = 95 GeV and (b) as a function of for 

minimal and maximal coupling to the Z boson. 



Assuming < mj> the main decay modes of t\ are ti —> ex® and t± — > The latter 

decay has practically 100% branching ratio if it is kinematically allowed. As \t further decays 
into Xi or Xi <?<?' ^e signature is two acoplanar b jets accompanied by two charged leptons + 
large missing energy (If), or single lepton + jets + 1$, or jets + If. Here the b tagging technique 
can be used to extract the signal. However, in this case the Xi w iU most likely be discovered 
first. On the other hand, if < m~+ +m b the flavour changing decay t\ ex? has practically 
100% branching ratio. (The decay t\ — > bff'xi, proceeding via a virtual chargino, is negligible 
35f.) The signature is then two acoplanar jets + If. Quite generally, the invisible energy is 



larger for t\ —* cxi than in case of t\ — > bxt- If, however, m^++mb > > mj+^+m^+m^) 

the decays t\ — > bvt^ or t\ — > b£ + i), proceeding via a virtual xt i can compete with the decay 
into cxi 0, 0, iH- In this case the signature is 2b + 2£ + If (or b + £ + + jet + If or 



jets + $). 

The main decay modes of &i are &i — > fo^i an d &i — * &X2> the second decay being possible in 
the parameter region approximately given by M < m~ bi — or < — m&. For the 

channel the signature is two acoplanar b jets + missing energy If. If the b± decays into bxl the 
b jets are acompanied by additional jets and/or leptons from x\ ~^ XiQQ and/or x% ~^ Xi$- b 
tagging will help to enhance the signal. 

If the lifetime of the squark is longer than the typical hadronization time of O(10~ 23 s), i. e. 
T < 0.2 GeV, it hadronizes first into a colourless (qtf) or (q x qq) bound state before decaying 
3T|. This is generally expected in case of t\ — > cxi and t\ — > bu£ + , bi + v since these decays 



involve the electroweak coupling twice ||35|| , but can also occur for i\ — > bxt an d &i ~ ^ t>Xi 
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depending on the nature of \i an d X?, and the squark mixing angles p8|, see also Chapter 3 



2.1.2 Linear Colliders with y/s = 0.5 - 2 TeV 

A future e + e~ Linear Collider will in many aspects be complementary to the CERN LHC. Its 



physics capabilities have been studied in ^4], [55], [56], [57| , see also |23], ^4]]. The virtues of an 
e + e~ Linear Collider are: 

• a very clean environment, 

• flexible centre-of-mass energies, and 

• high polarization of the e~ and possibly also of the e + beams. 

One can thus tune an e + e~ Linear Collider for many different purposes. For instance, one 
can optimize -y/s for specific production processes. Moreover, one can use beam polarization to 
enhance signals and to suppress backgrounds. One can also make complementary measurements 
by using both polarization states of e~ and e + . This may be of great advantage for testing the 
TeV range in a conclusive form. 

There are various Linear Collider projects under study in Europe, Japan, and the USA |58| . A 
majer breakthrough was achieved recently when it was realized that an integrated luminosity 
of £ = 300 fb" 1 at = 500 GeV and of C = 500 fb~ x at = 800 GeV can be reached with 
the TESLA design (assuming a Snowmass year of 10 7 s for running). 



In this section, we present the cross sections of stop and sbottom pair production at an e + e~ 
Linear Collider for unpolarized beams, polarized e~ beams, and polarized e + and e~ beams. 
For squark production in e + e~ annihilation with unpolarized beams and polarized e~ beams 
see also |39|, |60l . 



The dependence of the e + e — > t\t\ cross section is shown in Fig. |2l|a for = 
180 GeV, cos6^ = 0.7 and unpolarized beams. The effects of SUSY-QCD corrections from 



gluon and gluino exchange and of initial state radiation are demonstrated in Fig. |2]3|b. Note 
that at high energies the gluino exchange contribution has the opposite sign of the gluon 
exchange contribution, and the absolute values are increasing with y/s. The effect due to initial 
state radiation turns out to be of the order of 10%. The sum of all corrections can well exceed 
10%. The effects are similar for e + e~ — > b]b\. 
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Figure 2.3: (a) Total cross section and (b) radiative corrections relativ to the tree-level cross 
section for e + e~ — > t\ti as a function of \fs for = 180 GeV and cos# t ~ = 0.7 (m ( - 2 = 300 GeV, 
m- g = 300 GeV) 



Figure |2.4| a shows contour lines of the total e + e — > t\ii cross section in the — cos Of 



plane for y/s = 500 GeV and unpolarized beams. Analogously, Fig. [2.4| b shows contour lines 
of the total e + e~ — > b\b\ cross section in the m~ bi — cos6~ b plane. For the calculation of the 
SUSY-QCD radiative corrections we have used m 4 ~ 2 = mr = m g = 300 GeV. The &i&i cross 
section is about two to four times smaller than the tit± cross section. 

Analgous contour lines of ^2 and 62^2 production cross sections at y/s = 2 TeV are shown in 



Fig. [2.5| . Here we have used g = 300 GeV and nig = 700 GeV for the radiative corrections. 
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At the first stages of a Linear Collider pair production of the heavier mass eigenstate t 2 
(b 2 ) may be out of reach. However, the associated production of t\ and t 2 (&i and b 2 ) may be 
possible. Figure ^6]a shows the cross sections of t±t 2 production as contour lines in the m% - 

cos6i plane for rrif = 600 GeV and y/s = 1 TeV. In Fig. |2.6| a (e + e~ — > &1&1) at y/s = 1 TeV is 
shown as a function of cos 6~ b for (mr , m b2 ) = (300, 320) GeV and (400, 450) GeV. In both plots 
nig = 600 GeV. As one can see, the cross sections are large enough to be detected in a large 
region of the accessable parameter space. Here notice that replacing the stops by sbottoms 
or vice versa in Fig. |2.6| one gets essentially the same pictures. The only differences are due 
to SUSY-QCD corrections. Whether or not one can distinguish the two mass eigenstates of 
course depends on their decay properties. 



400 F 




0. 0.2 0.4 0.6 0.8 1. 0. 0.2 0.4 0.6 0.8 1. 

cos 61 cos 6~ b 



Figure 2.6: (a) Iso-cross section lines for t\t 2 production in the mj - cos 6* f - plane for y/s = 1 TeV, 
mi 2 = 600 GeV, and m g = 600 GeV. (b) Cross sections of e + e~ — > b\b 2 at y/s = 1 TeV as 
a function of cos^ for (m~ bl , mjj = (300, 320) GeV and (400, 450) GeV. Both plots are for 
unpolarized beams. 



Let us now turn to the effects of beam polarization. Figure ^7\a. shows the cos 61 dependence 
of the e + e~ — > t\ti cross section for left- and right-polarized as well as for unpolarized e~ beams 
for y/7s = 500 GeV and m h = 180 GeV (m t ~ 2 = m g = 300 GeV). For both left- and right- 
polarized e~ beams the cross sections depend strongly on the mixing angle. It is important to 
note that this dependence is opposite for left and right polarization. Therefore, experiments 
with polarized e~ beams will allow a more precise determination of the mass and the 
mixing angle 9$ from cross section measurements. This issue will be discussed in detail in the 
next section. Similar arguments hold for b\bi production with polarized e~ beams as shown in 
Fig. [2.7| b. However, in this case the dependence on the beam polarization is less pronounced 
for cos6^ < 0.5. The influence of beam polarization on t 2 t 2 and b 2 b 2 production at y/s = 
2 TeV is demonstrated in Fig. ^]8| for ^ = 700 GeV. Again, the cos 9g dependence is much 
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stronger for polarized than for unpolarized beams, however, the behavior is opposite to that 
of e + e~ — > q 1 q 1 . For the calculation of SUSY-QCD corrections we assumed g = 300 GeV 
and m-g = 700 GeV. 
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Figure 2.7: cos^ dependence of the cross section of (a) e + e — > titi and (b) e + e — > &1&1 
for left- and right-polarized as well as for unpolarized beams for = 500 GeV and 



m h ii = 180 GeV. 
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Figure 2. J 
left- and 
700 GeV. 



V. cos6g dependence of the cross section of (a) e + e — > £2^2 and (b) e + e — > 62^2 for 
right-polarized as well as for unpolarized e~ beams for = 2 TeV and m t - 2 1 2 = 
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If both beams are polarized these effects are even enhanced. This is shown in Fig. |2.9| where 
we plot the cross sections of t\ and b\ pair production as functions of cos 9g for yfs = 500 GeV, 
m ti h = -^0 GeV, and various beam polarizations: The full lines are for 90% polarized e + and 
e™ beams, (V_,V + ) = (—0.9,0.9) and (0.9,-0.9); the dashed lines are for 90% polarized e~ 
and 60% polarized e + beams, (V_,V + ) = (—0.9,0.6) and (0.9,-0.6); the dotted lines are for 
90% polarized e~ and unpolarized e + beams, (V_,V + ) = (—0.9,0) and (0.9,0). For the SUSY- 
QCD corrections we have taken m 4 ~ 2 = mg = m g = 600 GeV. Analogously, Fig. |2.10| shows 
the cross sections of &1&2 production at yfs = 500 GeV, for mg = 220 GeV, mg 2 = 240 GeV, 
m g = 600 GeV, 90% polarized e~ beams, and 90% and 60% polarized as well as unpolarized 
e + beams. The dependence on the degree of polarization can be seen in Fig. [2.11j : Here we 
plot er(e + e~ — > t\ti) and a(e + e~ — > &1&1) as function of P = V_ = —V + for y/s = 500 GeV, 



m 



I bi = 180 GeV, and the three cases q 1 = q L , q l = q R , and q ± ~ |(g L + q R ) 
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Figure 2.9: cos^ dependence of the cross section of (a) e + e~ — > titi and (b) e + e~ — > bibi for 
yfs = 500 GeV, mi g = 180 GeV, and 90% polarized e~ beams. The solid (dashed) lines are 
for \V+\ = 0.9 (0.6) and sign(7 :, + ) = — sign("P_). The dotted lines are for V + = 0. 
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Figure 2.10: cos^ dependence of the e + e~ — > &1&2 cross section for ^Js = 500 GeV, = 
230 GeV, m~ b2 = 260 GeV, and rrig = 600 GeV. The solid line is for unpolarized beams, the 
dashed lines are for \V-\ = 0.9 and V+ = 0, and the dashdotted lines are for \V-\ = \V+\ = 0.9 
and sign(P+) = -sign("P-). 
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2.1.3 Determination of Soft— Breaking Parameters — A Case Study 

In this section we estimate the experimental accuracies that can be obtained for stop and 
sbottom masses and mixing angles at an e + e~ Linear Collider. 

One possible way to determine and cos 9$ is using the y/s dependence of the unpolarized 
e + e~ — > t\t\ total cross section. Let us take mj = 180 GeV and cos% = 0.57 as reference 
point, and \fs = 400 GeV and yfs = 500 GeV as the two reference energies. Note that at 
| cos | ~ 0.57 the t\ti cross section has its minimum. Motivated by the Monte Carlo study of 
39[ |6T| and the recent big progress in the luminosity |59fl , we assume that the titi production 



cross sections can be measured with accuracies of ±5% at \/s = 400 GeV and ±2.5% and at 
yfi = 500 GeV. This leads to: 

du = 26.2 ± 1.3 fb at ^ = 400 GeV, (2.11) 
a v = 59.1 ± 1.5 fb at = 500 GeV. (2.12) 

For the SUSY-QCD corrections we have used m ( - 2 = 285 GeV and rrig = 500 GeV. Figure [2.12| a 
shows the corresponding error bands and the error ellipse in the m^-cos^ plane. As can be 
seen, the mass of t\ can be determined with good precision with this method. However, the 
precision on the mixing angle is rather poor. 

The polarization of the e~ beam offers the possibility of measuring the sfermion mixing 
angles with much higher accuracy. The cross sections of e + e~ — ► t\t\ for 90% left- and right- 
polarized e~ beams at the reference point m% x = 180 GeV and cos#f = 0.57 for yfs = 500 GeV 
are 

a_ 09 = 61.2 ±1.5 fb, (2.13) 
a +0 . 9 = 57.1 ± 1.4 fb, (2.14) 

where we have assumed that an accuracy of ±2.5% can be achieved. Figure |2.12| b shows 
the correponding error bands and the error ellipse in the m^-cos^- plane. The experimental 
accuracies obtained in this way for and cos 9f are: 

m h = 180 ± 1.65 GeV, (2.15) 
cos f = 0.57 ±0.012. (2.16) 

We next treat the sbottom system in an analogous way. If tan/3 is not too large we can 
neglect left -right mixing in the sbottom sector. From RGEs we expect b\ = hi and 62 = 
in this case. As reference point of the sbottom system we therefore take m~ bi = 200 GeV, 
mi 2 = 220 GeV, and cos6^ = 1. Based on |39| and we assume that the cross section 

for e + e~ — ► 6161 with 90% left-polarized e~ beams can be determined with an experimental 
accuracy of ±2.5%. For e + e~ — ► b 2 b 2 with 90% right-polarized e~ beams we take ±8% as the 
experimental accuracy. This leads to: 

a_ og (e + e- -> hh) = 82.3 ± 2.1 fb, (2.17) 

CT +o.9 ( e+e ~ hh) = 8.6±0.7fb. (2.18) 
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Figure 2.12: Error bands and error ellipses for the cross section of e + e~ — > t\ti in fb as a 
function of rrif and cos#f, (a) for unpolarized beams at y/s = 400 and 500 GeV and (b) for 
90% left- and right-polarized beams at \/s = 500 GeV. The dot corresponds to the reference 
point = 180 GeV and cosflf = 0.57. 



Here we have assumed that sbottom and stop production can be distinguished. This is possible 
in the parameter region where the sbottoms mainly decay into bx\- However, in the parameter 
region where the decay into * s possible, signatures similar to those of stop production could 
occur. The errors for the sbottom masses follow as: 

m 5l = 200 ±1.03 GeV, (2.19) 
m- b2 = 220 ± 1.93 GeV. (2.20) 

With these values for , cos 9$, mg , and m~ b2 we can use Eq. ( |1.51| ) and obtain the mass 
of the heavier stop t 2 if tan (3 is known from other experiments. Taking, for instance, tan (3 = 3 
leads to 

m h = 285.3 ± 3.3 GeV. (2.21) 

Confirming this value by t±t2 and/or £2^2 production at higher energies would be an independent 
test of the MSSM. 

Assuming that also // is known from an other experiment we are now able to calculate the 
underlying soft SUSY breaking parameters Mq, Mfj, M^, A t and Ab for the squarks of the 
third family, see Eqs. (1.46) to ( |1.50|) . Taking e.g., \i = 300 GeV (and tan/5 = 3) we obtain the 
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following values (A b = //tan/? = 900 GeV since cos^ = 1): 

Mq = 192.8 ±1.1 GeV, (2.22) 

Mjj = 136.9 ± 5.5 GeV, (2.23) 

Mq = 218.8 ± 1.9 GeV, (2.24) 

A t = -31.2 ±8.1 GeV. (2.25) 



While Mq, Mjy, and are obtained with good accuracy, the error on A t is quite large. This, 
however, depends on our specific choice of cos%, tan/3, and /t. For cos#£ = —0.57 (tan/3 = 3 
and fi = 300 GeV), for instance, we get m h = 180 ± 1.63 GeV, A t = 231.2 ± 8.1 GeV, and the 
other parameters as above. 

An alternative way to determine the squark masses and mixing angles is the kinematical 
reconstruction of the squark decays. This has been studied in e.g., p2"| , and for the 
case t\ — > bxt ■ 

As a last remark we note that the high precision that can be expected at a Linear Collider 
indicates that radiative corrections should be taken into account in the Monte Carlo studies. 
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2.2 Cross Sections for Muon Colliders 



The idea of a muon collider has already been brought up in the 1960 's and 70 's by Tinlot, 
Skrinsky, Neufer, et al. However, major technical challenges have to be met in order to collect, 
accelerate, and collide muons and antimuons, which decay with a lifetime of 2.2 /is. Not until 
1995, at the Sausalito workshop, was it realized that with modern ideas and technologies a muon 
collider may be feasable. The next step further was the Muon Collider Feasibility Study for the 
Snowmass 1996 workshop pS] . At present, the Muon Collider Collaboration carries out the 



R&D for a muon collider in the US. In Europe, a muon collider is discussed within the ECFA 
Prospective Study on Muon Colliders |67|] , see also [BI|. Specific studies are done for 0(100) 
GeV, 0.3 — 0.5 TeV, and multi-TeV machines with luminosities of 10 31 to 5 x 10 34 cm~ 2 s~ l . 

The advantages of a muon collider are the following: 



• The effective energy of a lepton collider is much larger than that of a hadron collider with 

the same center-of-mass energy. 

• In contrast to electrons, muons generate almost negligible synchotron radiation |. 

• The cross section for direct (s-channel) Higgs production in lepton-antilepton annihi- 
lation is proportional to mf. It is thus 40000 times larger at a than at an e + e~ 
collider. 

• Because of the lack of bremsstrahlung (and synchotron radiation) energy spreads as small 
as 0.003% are expected. By measuring g-2 of the muon it should be possible to determine 
the absolute energy to an even higher accuracy. 



A muon collider can thus be circular and much smaller than e + e~ or hadron colliders of compa- 
rable effective energies. With its expected excellent energy and mass resolution a muon collider 
offers extremely precise measurements. Moreover, it allows for resonant Higgs production; in 
particular it may be possible to study the properties of relatively heavy H° and A , which can 
hardly be done at any other collider. 

The technical challenges still to be met include (i) a liquid metal target for pion production, (ii) 
a new cooling technique (ionization cooling?) for the muons, and (iii) the problem of neutrino 
radiation j70fl . For a more detailed introduction, see e.g. ffl], [f2]| . Information about ongoing 



work can be found in |6f 

Squark pair production in fi + fi~ annihilation proceeds via the exchange of a photon, a Z 
boson, or a neutral Higgs boson. The corresponding Feynman diagrams are shown in Fig. p. 13 . 



Provided that m H > nig- + rriq. (H = H°, A°)Q the process fi + fi — > q^g- thus offers the 
possibility to study the squark-squark-Higgs couplings at the Higgs boson resonances. 



3 The energy loss per revolution in a circular machine is proportional to m 4 



/■ 

4 As the MSSM predicts m h o < 130 GeV squark production at the h° resonance is already excluded. 
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Figure 2.13: Feynman diagrams for the process ji + ji — > q^j, (a) i = j, (b) i ^ j. 



by 



For unpolarized beams the differential cross section at tree-level (to order m^/y/s) is given 



d a 37ra 2 K,j 



dcosi? 4s 2 



Til 2 — Til 2 tc 

^ T vv sin 2 $ + T HH + a Tv \ + *L cos 



VH ' J-VH' 



[2.26) 



with the kinematic function = [(s — m~. — m~ ) 2 — 4m~ i m~ J 1 / 2 and i? the scattering angle of 

denotes the 7 and Z exchange contributions (which are the same as for e + e~ — > q^j), T HH 
denotes the contributions from Higgs boson exchange, Ty H comes from the Z- A interference, 
and Ty H from the (7, Z) - (h°, H°) interference: 



Tyy 



zb W L W 



16 SjyCjy 



(2.27) 



2e 4 



(Gi)ij sin « ^ - (Cl) i7 - cos ad H + sin /3 d A 



' 2.7 



(2.28) 



rpa 



m ^K sin Pc^Gl)^ 



V2 



e2 c 2 .-,2 



(2.29) 



rpb 



2\[2 eq-mphpSij 



(GDtj sin a - (Gf)^- cos a 9fJ(d#) 



myVnK Cjj s 
s w c w 



(Gf )„ sin a »(d^d fc ) - (Gf )„ cos a K(<£d H )l (2.30) 
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with Vn = —1+4 sin 2 9w and a M = — 1 the vector and axial vector couplings of the muon to the 



Z boson; c^- is the squark coupling to the Z, Eq. ( |1.95| ), G 9 ^ 23 j that to {h°, H°, A } as given 
in Section 11.6.4! FL and 



d 



x 



[s - m x ) + iT x m x 



-i 



(2.31} 



The total cross section is given by 



a (/i + /U~ 



3ira k 



2s 2 



2^ 
3 s 2 



HH 



2 2 

m~ — m~ 

' VH 



(2.32) 



Notice that the 7 exchange diagram does only contribute to q t q { production. On the other 
hand, since {G\)n = 0, A exchange only contributes to q x q 2 production. Notice moreover that 
the (7, Z) - (h°, H°) interference in ( |2.26| ) is proportional to cost?, giving rise to a forward- 
backward asymmetry. Being proportional to m M this asymmetry is, however, numerically very 
small, typically of the order of 10~ 4 or less. 

Away from the Higgs boson resonances one has a(fi + fi~ — > g^g.) ~ a(e + e~ — > q^qj). However, 
for \fs ~ m j^o ^0 the cross section at a muon collider is largely enhanced compared to that of an 
e + e~ collider due to resonant H° and/or A exchange. While the pure gauge boson contribution 
shows a sin 2 $ angular dependence, the pure Higgs boson contribution does not depend on the 
scattering angle. This can be used to disentangle these two contributions. 



However, since the muons originate from pion decays, they are naturally polarized. The 
differential cross section for a fi~ beam with polarization V- and a u, + beam with polarization 
V+ is: 



da P 
d cos^ 



3tTQ> Ki 

4s 2 



s 2 



vv 



sin 2 $ + T HH 



+ % 



VH 



(2.33) 



with 



el 8. 



q u l3 



'1 - vn 



+ 



c 2 s 2 

16 Sy^Cyy 



(1 - V-V+) - a M (V- - V + )\ K(d z 



[vl + aj) (1 - VJP+) - 2 v a„ (V- - V + 



\d z \\ 



(2.34) 



2e 4 



(Gi)ii sin a d h ~ ( G l)ij cos " d H " + (GDa sin P d A 



A + vn 



+ 2 sin (Gfyij (Gf k- sin a $t(d* A d h ) - (Gf cos a &{d* A d H ) (V- + V. 



(2.35) 



5 You may have noticed that (C7g) 12 = —{G\) 2 i- Fortunately, this difference in sign is "repared" by the term 
in front of Ty H and thus cr(/i + /i _ — > Q1Q2) = cr (/ i+ M _ — * Q2Q1) as it should be in case of unpolarized 
beams. 
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and 
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(2.36) 



V± = (— 1, 0, 1) for (left-polarized, unpolarized, right-polarized) /x ± beams. The total cross 
section is given by: 



o>(//V -> 



37ra 



2 lyv + Ai 



(2.37) 
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(2.38) 



where = /cy/s and 



Z,/ v^s^c 2 



m=. — m 



sin/? (G1V. -£(?-- V. 



M(d* z d A ) 



+ y 



(Gf) tf sin a &(d* z d h ) - (G|) i7 - cos a 3J(d^d H ) (P- + 7> 



(2.39) 



Note the additional if ) interference in case of q x q 2 production. This interference term 

changes sign when one goes to the CP conjugate state. Analogous terms occur for the Z- 
(h°, H°, A ) interferences. Hence, cx(/i + /i~ — > qiq 2 ) ^ cr(/i + /i~ — > q 2 q x ) if V_ ^ —V + ! 
Note moreover, that /i^/i^ and p^fiR combinations only couple to Higgs bosons while [i^nl 
and fJ-t^R combinations only couple to gauge bosons. Therefore, for Higgs- and related physics 
it is highly preferable to collide fi + and yT of same helizities. It is expected that a polarization 
of 28% for both beams — and maybe more with loss of luminosity — can be achieved (see 
e.g. @). 



In Fig. 2.14 we show the dependence of the total tii\ production cross section for = 
180 GeV, cos6>j = —0.55, and tua = 450 GeV. The three plots are for unpolarized (P = V_ = 
V + = 0), 30% left-polarized (P = -0.3), and 60% left-polarized (P = -0.6) fj+ and yr 
beams, respectively. In addition to the cross sections at a yU + /i~ collider (full lines), the dashed 
lines show the cross sections in e + e~. For calculating the properties of the Higgs bosons we 
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have taken m ( - 2 = 260 GeV, m bi = 175 GeV, m. 



& 2 



0.9, M L = 170 GeV, 



150 GeV, A T = 300 GeV, M = 140 GeV, y 



parameters we get m H o = 454 GeV, sin a = —0.35, T 



195 GeV, cos6- b 
300 GeV, and tan/5 = 3. With this set of 
H o = 5.4 GeV, and I^o = 7.3 GeV. As can 
be seen, a clear peak occurs at y/s = m H o . The total cross section at the resonance varies only 
little with the polarization. However, the relative importance of the Higgs exchange strongly 
increases with increasing polarization. 

Analogously, Fig. [2.15| shows y/s dependence of <j(/i + /i~ — > tii 2 ) and a(ji + n~ — > i 2 £i) for the 
parameters of Fig. [2.14j . In this case also A exchange contributes. While <7(tit 2 ) = cr(^2^i) for 



unpolarized beams, there is a difference between these cross sections for P / 0. This difference 
strongly increases with increasing polarization. Notice, that while for t x t 2 production the A 
and H° resonances fully overlap this is not the case for t 2 t\ production. 

This can be seen more clearly in Fig. J2.16] which zooms on the resonances of Figs. |2.14| and 
|2.15| : <j(t\ti) peaks at y/s = m#o, a (tit 2) has its maximum at y/s ~ itla, and <r(i 2 ti) shows two 
humps at y/s ~ rriA y H - With the excellent energy resolution of a muon collider it may thus 
be possible to determine the properties of heavy neutral SUSY Higgs bosons. In general, it is 
expected that this will be difficult at hadron and e + e~ colliders. 



1.12 Mq, A t 
300 GeV, tan/3 



A b , M L = 170 GeV, M E = 150 
3, and itla = 450 GeV. In 
454 GeV. The dashed lines 



The cos Qi dependence of a(y + fi — > t\ti) is shown in Fig. |2.17| for 30% left-polarized ji beams 
and m h = 180 GeV, m h = 260 GeV, M b 
GeV, A T = 300 GeV, M = 140 GeV, y. 
Fig. [TIT] a y/s = m A = 450 GeV and in Fig. pT7| b yfs = m H o 
show the Higgs boson contributions, the dotted lines the (7, Z) contributions, and the full lines 
the total cross sections. While the vector boson contributions depend only on cos 2 Qi the Higgs 
boson contributions induce a distinct dependence also on the sign of the stop mixing angle. 
Therefore, the cos^ dependence of the t\ii production cross section at a ji + n~ collider with 
yfs ~ rriA is very different from that at an e + e~ collider. 

The same is true for tit 2 production. This can be seen in Fig. |2.18| where we plot a(n + fi~ — > tit 2 ) 
and a(y + fi~ — > t 2 ^i) as a function of cos 6^ for the parameters of Fig. |2.17| . Here note, that the 
leading term of the t^H coupling is hf(y cos a — A t sin a) cos 2% while the t^A coupling is 



h t (n sin (3 — At cos (3). Since /1 is constant in Fig. |2.18| , A t varies with cos 9 



An example for sbottom production is shown in Fig. 2.19| where we plot the yfs dependence 
of the cross sections of 6161, &i& 2 + b 2 bi, and 6 2 6 2 production for 30% left-polarized y^ beams. 
The input parameters are: m 



mi 2 = 256 GeV, cos^- 



6j 

-0.66, M 



157 GeV, m- b2 
160 GeV, 



188 GeV, cos9 b = 0.78, m h = 197 GeV, 



GeV, y, = 300 GeV, m A = 380 GeV, and tan (3 = A. 



= 155 GeV, A T = 100 GeV, M = 140 
In this case we get m H o = 383 GeV, 



sin a = —0.28, T H0 = 1 GeV, and r A0 = 1-9 GeV. Note, that for bib 2 production (Fig. |2.19| b) 



two distinct peaks occur at the A and H° resonances! However, for the parameters of Fig. [2.19| 
it might be difficult to distinguish bi and 6 2 because their decay properties are very similar. We 
therefore show in Fig. [2.20| a the total sbottom production cross section, Y^ij cr(y + y~ — > hbj), 
for the parameters of Fig. [2.19| . In addition to P = —0.3 (full line) we also show the case 
P = —0.6 (dashed line). As can be seen, the two peaks are still separated. Figure p.20| b 
compares the 6i6 2 and & 2 &i production cross sections for P = —0.6 and the other parameters 
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as in Fig. [2.20| a. For P = —0.3 the difference between a{pib 2 ) and a{b 2 bi) is less pronounced. 
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Figure 2.14: Cross sections of fi + fi — > titi for unpolarized, 30% and 60% left-polarized fj^ 
beams as a function of y/s, for = 180 GeV, m t - 2 = 260 GeV, cos^- = —0.55, m~ b = 175 
GeV, m- b2 = 195 GeV, cos0 s = 0.9, M L = 170 GeV, Mg = 150 GeV, A T = 300 GeV, M = 140 
GeV, [i = 300 GeV, tan/5 = 3, and m A = 450 GeV; P = V_ = V + . 
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Figure 2.15: Cross sections of ji + n — > t\t 2 and n + \i — > t 2 t\ for unpolarized, 30% and 60% 
left-polarized beams as a function of y/s, for = 180 GeV, mj 2 = 260 GeV, cos#f = —0.55, 
m - bi = 175 GeV, m b2 = 195 GeV, cosflg = 0.9, M L = 170 GeV, M s = 150 GeV, A T = 300 GeV, 
M= 140 GeV, n = 300 GeV, tan/3 = 3, and m A = 450 GeV; P = V_ = V + . 
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Figure 2.16: Cross sections of fi + fi — > ttfj for 30% left-polarized fj^ beams (V_ = V + = —0.3) 
as a function of y/s, for m~ tx = 180 GeV, mi 2 = 260 GeV, cos#£ = —0.55, m~ bi = 175 GeV, 
m- b2 = 195 GeV, cosflg = 0.9, M L = 170 GeV, Mg = 150 GeV, A T = 300 GeV, M = 140 GeV, 
H = 300 GeV, tan/5 = 3, and m A = 450 GeV. 
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Figure 2.17: Cross section of — > titi for 30% left-polarized beams as a function of cos# t ~, 
for m h = 180 GeV, m h = 260 GeV, Mf, = 1.12 Mq, A t = A 6 , M z = 170 GeV, = 150 
GeV, A T = 300 GeV, M = 140 GeV, n = 300 GeV, tan/5 = 3, and m A = 450 GeV; (a) 
y/s = 450 GeV, (b) y/s = 454 GeV. The dashed lines show the Higgs boson contributions, the 
dotted lines are the (7, Z) contributions, and the full lines are the total cross sections. 
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Figure 2.18: Cross sections of fi + fi — > tfa and fi + fi — > £2^1 for 30% polarized beams as a 
function of cos% for m~ tl = 180 GeV, m h = 260 GeV, M D = 1.12 Mq, A = A b , M z = 170 
GeV, = 150 GeV, A T = 300 GeV, M = 140 GeV, fx = 300 GeV, tan/5 = 3, and m A = 450 
GeV; in (a) = 450 GeV and in (b) = 454 GeV. 
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for 30% left-polarized beams (V_ = V + 
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Figure 2.20: (a) Total sbottom production cross section and (b) cross sections of — > b^b 2 
and — > b 2 b\ as a function of y/s, for = 157 GeV, m~ b2 = 188 GeV, cos9~ b = 0.78, 
m h = 197 GeV, m h = 256 GeV, cos# t ~ = -0.66, M~ L = 160 GeV, M% = 155 GeV, A T = 100 
GeV, M = 140 GeV, fx = 300 GeV, m A = 380 GeV, and tan [3 = 4. 



Chapter 3 

SUSY-QCD Corrections to Squark 
Decays 

3.1 Introduction 

We now turn to the stop and sbottom 2-body decays. As the squarks of the 1st and 2nd 
generation, stops and sbottoms can decay into quarks plus charginos, neutralinos, or gluinos. 
In addition to these "coventional" decays into fermions, stops and sbottoms may decay into 
bosons, i.e. into a lighter squark plus a gauge or Higgs boson. These decays are possible in 
the parameter space where large Yukawa couplings and q^-q^ mixing generate a large mass 
splitting between q x and q 2 and/or q t and qL Moreover, Yukawa couplings and q L -q R mixing 
have a strong influence on the branching ratios of the various decays modes. Thus, stops and 
sbottoms can have quite complex decay patterns. 

In summary, there are the decays (i, j = 1, 2; k = 1 . . . 4): 

bxl,txj, (3-1) 

bg, (3.2) 

ij + (W-,H-), (3.3) 

h + (Z°, h°, H°, A ). (3.4) 



U -> txl, bXj, bi 

U -> tg, bi 

U -> % + (W + , H + ), b t 

t 2 -> h + (Z°, h°, H° } A ), b 2 



All these decays were first discussed at tree level in J74 



in 



They were studied in more detail 
and b\ the decays into fermions of 



C| |6T], ^9], |75| . For the lighter mass eigenstates t\ 
Eqs. (3.1) and (3.2) are in general the most important ones []. If the strong decays into gluinos 
are kinematically allowed they play an important role. For the heavier £ 2 and 62 it turned out 



1 If, however, all 2-body decay modes are kinematically forbidden, loop-induced decays such as <i — ► 
and 3-body decays become important. The higher order decays of t\ have been studied in [^5[ Ell ^3| . 



cx°i 
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the bosonic decays of Eqs. (3.3) and (3.4) can be dominant in a wide range of the parameter 
space; see [[fjj for a detailed analysis. Moreover, also the decays b\ — > t\W~ and b\ — > t\H~ 
can be important if the t-h mass splitting is large enough |76| . 



The 1-loop SUSY-QCD corrections to the decays of Eqs. (3.1) to (3.4) have all been calcu- 
lated within the last two years (1996-98). For the decays into charginos and neutralinos this 
was done in Jf7|, [78|, [79! ' f° r the decays into W ± and Z bosons in |80f , for the decays into Higgs 
bosons in fl8T] , |82|| , and for the decays into gluinos in [^, 75]. (Here note that the decays into 
photon or gluon, which are absent at tree level, are not induced by these corrections either.) 

In this chapter we study the 0(a s ) SUSY-QCD corrections to the decays of Eqs. (3.1), (3.3), 
and (3.4). In Sect. |3.2| we first discuss some general aspects of these calculations. Notations 



and conventions are clarified in Sect. |3.2.3| ; the formulae for the self-energies of squarks and 
squarks are given in Sect. |3~B| . In Sect. we discuss the subtleties that have to be taken into 
account for calculating the 0(a s ) SUSY-QCD corrections in the on-shell scheme. The specific 
squark decays are then treated in Sects. |3.5| to |3^. For each process, we give the complete 
formulae for the 0(a s ) SUSY-QCD corrected decay widths and perform a detailed numerical 
analysis. Finally, the SUSY-QCD corrected branching ratios are discussed in Sect. [Tl8[ 



3.2 General Aspects of SUSY-QCD Corrections 

For the corrections at 0(a s ) we have to consider 1-loop diagrams with gluons, gluinos, and 
squarks in the loops. The integration momenta in the loops run from zero to infinity, and some 
of the diagrams diverge. We therefore have to apply appropriate renormalization procedures 
to cosistently isolate and remove the infinities from the measurable quantities. The general 
approach is as follows: 



3.2.1 Renormalization 

We start with the bare Lagrangian, which consists of bare parameters (bare masses and cou- 
plings) and bare fields, and write it in terms of the renormalized quantities. The bare mass 
m and bare coupling Co, for instance, are replaced by the renormalized mass m and coupling 
c plus the associated counterterms: 

m = m + 5m, Co = c + 5c. (3.5) 

The bare fields are given by the renormalized ones multiplied by a wave-function renormaliza- 
tion factor e.g., 

O = Z 1 ^ (3.6) 

(where Za is in general a matrix which mixes fields with the same quantum numbers into each 
other). With these replacements we can seperate the bare Lagrangian into the renormalized 
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one and a part which contains the counterterms: 



o — 



C + 5C. 



(3.7) 



Here we note that £ and 5C have the same structure. 5C gives rise to new (effective) vertices 
in each order of pertubation theory. Finally we apply appropriate renormalization conditions 
which determine the counterterms and the physical meaning of the renormalized parameters. 

For our calculations we use the on-shell renomalization scheme. Within this scheme the 
renormalized mass is defined as the pole mass, i. e. the experimentally measured mass 0. The 
on-shell renormalization conditions are: (i) the renormalized mass is the real part of the pole 
of the propagator and (ii) the real part of the residue of the pole is unity. 

Technically speaking, at first order we have to add 1-loop vertex and wave-function corrections 
to the tree-level diagram as well as the counterterms for the couplings. For the latter we have 
to renormalize all parameters that enter the couplings. The resulting 1-loop corrected decay 
width is then utraviolet (UV) finite. In order to cancel also the infrared (IR) divergence we 
add the emisson of real (hard and soft) gluons. 

The 0(a s ) SUSY-QCD corrected squark decay width can thus be decomposed in the following 
way: 



Here T is the decay width at tree level and ST^ V ' W ' C ^ are the virtual corrections. The superscript 
v denotes vertex corrections, w wave function corrections, and c the shift from bare to on-shell 
couplings. ST rea i is the correction due to real gluon emission. 



The divergencies which arise in calculating the loop integrals must be regulated. In the SM the 
standard technique is dimensioal regularization (DREG) which respects the gauge symmetry 
of the Lagrangian (and thus also preserves the Ward identities). In this scheme, spacetime is 
continued to D = 4 — e dimensions. For D < 4 the divergencies then appear as simple poles in 
e. However, in DREG also the dimensionality of the fields is continued to D dimensions. As a 
consequence, the index of a gauge field runs from to D — 1, and also the Dirac algebra is 
/^-dimensional e.g., 7^7 M = D. 

This causes problems in supersymmetry: SUSY requires that the numbers of bosonic and 
fermionic degrees of freedom be equal in each supermultiplet. Continuing a four-dimensional 
multiplet to D = 4 — e dimensions spoils this equality, introducing a mismatch of gauge boson 
and gaugino degrees of freedom. It also spoils the SUSY- Ward identities. 

2 This has to be distinguished from the minimal substraction method where the counterterms are 'purely 
infinte', i.e. they solely contain the divergent parts. In contrast to the on-shell scheme these counterterms are 
scale dependent. 



r = r° + 5T {v) + (5r (u,) + 5r (c) + sr. 



real • 




3.2.2 Regularization: DREG vs. DRED 
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A solution known as dimensional reduction (DRED) was proposed by Siegel, Capper, Jones, 
and van Nieuwenhuizen [[16], ITfl . In DRED spacetime is continued to D = 4 — e dimensions 
while the fields are not affected. Thus the index of the gauge fields runs from to 3, and 
the Dirac algebra is performed in four dimensions. This method nicely preserves both, gauge 
symmetry and SUSY — at least up to 2-loop order fl8[ . 



3.2.3 Notations, Conventions, etc. 



In this thesis we work in the on-shell renormalization scheme using dimensional reduction. For 



the explicit calculation of the loop integrals we use Passarino-Veltman functions |55] in the 
notation and convention of A. Denner 



86 



The conceptual and calculational details of renormalization in supersymmetry have been thor- 
oughly discussed in Refs. [87], |4"2f . We therefore here just present the final results of our 
calculations. Analytical expressions of Passarino-Veltman integrals, generic diagrams, as well 
as details on how to calculate individual graphs, are given in H42 |. 



For the numerical analysis we use on-shell squark masses 2 and mixing angles 9q (0 < 
9q < 7r) as input parameters. Moreover, we take m t = 175 GeV, mb = 5 GeV, mz = 91.2 
GeV, sin 2 9w = 0.23, a{mz) = 1/129, and a s {mz) = 0.12. For the running of a s we use 
Oi s {Q 2 ) = 127r/[(33 — 2nj) ln(Q 2 /A^)] with nj the number of flavors. We take a s = a s (mqj 
for the q { decay except for (i) the renormalization of Mq for which we take a s = a s (Mg(t)), see 
Sect. |3.4.2| , and (ii) the calculation of the gluino mass, for which we take m- g = a s {m^)jaiM 
with five steps of iteration. For the radiative corrections to the h° and H° masses and their 
mixing angle a (— f<«< : |by convention) we use the formulae of [26]; for those to we 
follow Q. 



In order to respect the experimental mass bounds from LEP2 J49[ and Tevatron [[88| we impose 



> 85 GeV, m h o > 90 GeV, m-+ > 95 GeV, and rrig > 300 GeV. Moreover, we require 

5p (t-b) < 0.0012 [^] from electroweak precision measurements using the one-loop formulae 
of and A 2 t < 3 (M|(t) + M? + mfj, A\ < 3 (M? (6) + M? + m 2 Hl ) with m 2 Ha = {m\ + 



m 



cos 2 P — I m z and 



in 



Hi 



[m. 



+ m z ) sin 2 /3 — 2 m l [EH f rom tree-level vacuum stability. 



3.3 Self Energies to 0(a 8 ) 



3.3.1 Squark Self Energy 



The squark self energy to 0(a s ) gets contributions from gluon, gluino, and squark loops. The 



relevant Feynman digrams are shown in Fig. 3.1. 



3 Notice that Eq, p7[ have the opposite sign convention for the parameter ji. 
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Figure 3.1: Squark self-energy diagrams to 0(a s ). 



91,2 



1j 



The gluon-squark loop leads to 



-(g), 



2a 



m- 



) = S ~mi {2B o (mL0, 



m 



(3.9) 



The gluon loop due to the qqgg interaction gives no contribution: It is proportional to A (\ 2 ), 
where A is the IR gluon mass, and vanishes for A — > 0. The contribution due to the gluino-quark 
loop is 



sf'K) =~~ iMm 2 q ) +m\ i B 1 (m%mlm\ 



4a 



3 7T 



+ \m\ + (-)* m~ g m q sin 26 q] S (m? ,m?,mj)}, 



(3.10) 



Eg>(<)=E£'(mS,) = 



(5)^2 



4 a, 



-m, 



>m, cos29 q B (ml,m 2 g, m 2 ), 



and the squark bubble leads to (i 7^ i') 



4* Vl) = ^ {cos 2 2^A (m|) + sin 2 2^A (m|J}, 



3tt 



a 



6/T 



Note, thatS®(m? i ) = S®(<). 



<5m? the shift from the bare to the on-shell squark mass, is given by: 



EW(mJ,) + ESf>(m?) + E®(m|,) 



(3.11) 

(3.12) 
(3.13) 



(3.14) 



The squark wave-function renormalization constants Z n i(qA are: 

*{E®'*>(mS,)} 



(9,9) 



m|, - m| 



(3.15) 
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with t u (m 2 ) = dE u (p 2 )/dp 2 \p2 
2a s 



3vr 



S (m? , 0, mfj + 2m? S (m? , A 2 , m? 



9i ' ' 



(3.16) 



2a, 
" 3tt 



d/2 2 2\ , / 2 2 2\ n / 2 2 2\ 

B o{m s .,m s , m q ) + (m- -m ? - m^olm^-^m-, m ? ) 
-2m q mg(-l) 1 sin 2^S (m|, m|, m 2 



(3.17) 



The four-squark intraction does not contribute to 5Zu because £® = 0. On the other hand, 
the off-diagonal squark wave-function renormalization constant 5Zu/ gets no contribution from 
gluon exchange diagrams because they do not mix q x and q 2 - Bare and renormalized squark 
fields (g 4 ° and are related by 



# = (1 + ft + , 

$* = (l + §<^)9* + ^$ 



(3.18) 
(3.19) 



3.3.2 Quark Self Energy 



The self energy of a quark gets contributions from the gluon and gluino loops shown in Fig. |3~2 



9 





Figure 3.2: Squark self energy diagrams to 0(a s ). 

5m q , the shift from the bare to the pole mass of the quark q, thus has two contributions: The 
gluon exchange contribution is 

2a s 



6m[ g) 



3 7T 



m„ 



B (mi, 0, m 2 ) - B x {m% 0, mi) - r/2 



and the gluino contribution is 



(3.20) 



Smf 



3^ 



m n 



D/2 2 2\ |D ^2 2 2 \ 

B x (m, , m~ g , m 5l ) + Bi {m q , m § , m & ) 



sin 26*g 



B (m 2 q ,m 2 § ,m 2 qi ) - B (m 2 q ,m 2 g,m 2 q2 )] } . 



(3.21; 



For the quark wave-function renormalization constants due to gluon exchange one gets: 



9 9 3 7T 



S + 2m 2 ,{B -B x ) -r/2 



(3.22) 
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with B m = B m (m q , A 2 , m 2 ), B m = B m (m q , A 2 , m 2 ); those due to gluino exchange are: 
5z l a 2«.f cog2 Q _ B , + gin2 g _ B 2 + m 2 r^i + jja + !^ sin ( £i _ fi* 

O 7T L 



(3.23) 



2 a, 



JZf 0) = ^ j sin 2 ^ + cos 2 ff B? 

H 6 IT 



2 i 2 

+ m ql 



Bl + Bf + ^sin20 9 ~(5 o 1 - J B ( 



where 5* 



-ELfm 2 mi, m|.) and 



quark field g° and the renormalized one g is 
g° = {l + \8Z L q P L + \8ZfP R )q, 
<f = q(l + \8Z L q P R + \8ZfP L ). 



B m (m 2 q , m\, m| 



m 



(3.24) 



The relation between the bare 

(3.25) 
(3.26) 



The parameter r in ( |3.20| ) and ( 3.22|) exhibits the dependence on the regularization scheme: 
r = in DRED while r = 1 in DREG (remember that only DRED preserves SUSY!). Here we 
note that in our calculations r does not cancel in the final results. Therefore, there are also 
numerical differences between the two regularization schemes. 



3.4 Subtleties 



All in all we have to consider O(10) loop diagrams per process. This may not sound much; 
however, there are some subtleties which have to be taken into account: 



3.4.1 Renormalization of the Squark Mixing Angle 



The stop and sbottom couplings depend on the respective mixing angles which therefore must 
be renormalized. The first proper renormalization prescription for the squark mixing angle 
was given in |^7|. There the counterterm 86 q was fixed in the process e + e~ — > ^g- such that 
it cancels the off-diagonal part of the squark wave-function corrections (diagrams h and j of 



Fig. [2.1| with k ^ i). The idea is the following: The variation of the squark-squark-Z coupling 



Cij gives 



Sc ij 



2 c, 



12 

'22 ~~ C H 



C 22 ~ C H 



-2c 



89 s 



(3.27) 



12 



The requirement that the SUSY-QCD corrections to e + e — > q^ be UV finite leads to: 

5c n = 2 Cvi 66 9 = -2 6Z 21 cu, (3.28) 

5 C12 = (c 22 - c u ) 58g = -8Z 21 c 22 - 8Z 12 cu, (3.29) 

8 c 22 = -2 cu 50g = -2 8Z 12 c 12 , (3.30) 
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where the SZ^ are the squark wave-funktion renormalization constants, see Eq. (|3.15| ). Using 
Eq. (3.29), i.e. the process e + e" — > g x g 2 , to fix the squark mixing angle we get (here and in the 
folowing £ 12 = 9ft(£ 12 )): 



SB* 



^(^J c 22 - S 12 (m|)c u 



(c 22 - cn)^ -m|; 



(3.31; 



86 q thus gets contributions from gluon and gluino exchanges, 56 q = 56 'i + • Explicetly: 



(5) 



ft. 



6tc m-j 



sin 46^ 



A (m?)-A (m?) 



5# = 



m- g m q 



3 iliCm? 



D/2 2 2\ d/2 2 2\ 

B Q (m q2 ,m~ g , m q )c u - B (m qi ,m- g , m q ) c 22 



(3.32) 
(3.33) 



This scheme was also applied for SUSY-QCD corrections to squark decays in |77| , pO 
well as to Higgs decays into squarks in W\ - We will use this scheme in what follows. 



as 



There are also other possibilities of defining the on-shell squark mixing angle. In |78| , for 
instance, 59 q was fixed such that the renormalized self energy of the squarks remains diagonal 
on the q 2 mass shell (i.e. using Eq. (3.30)). This leads to the condition 



so, m) 



S12K) 



(3.34) 



An analogous condition on the q 1 mass shell was applied in |8T]. In |79[ a scale-dependend 
definition was used: 



, 9 -(Q 2 )([79])- -i2 



Si 2 (Q 2 



(3.35) 



A process independent "democratic" approach is to take the arithetic mean of Eqs. (3.28) and 
(3.30), as done in pifl . 50^ is then given by 



^([84]) 



S 12 (ir4) + E 12 (m|; 
2(m 2 qi -m 2 q2 ) 



(3.36) 



This is, by the way, the only definition that works for processes with charginos in the loops 
(Yukawa correctins etc.). 



The differences between the various schemes are ultraviolet finite. In Fig. |3.3| we compare 59^ 
of @ 0, to that of our scheme g7J. Considering squark decays, 5^-(Q 2 ) (|7|) = 
or 59 q (Q 2 ) (J75)) = (0D- The numerical differences between the various schemes are very 
small, typically well below 1%. 
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Figure 3.3: Differences in 5cos# t ~ between [78, BlL 84[ and our scheme [37[] as a function of 



cos^-, for m h = 250 GeV, m t ~ 2 = 600 GeV, and m- g = 600 GeV. 
3.4.2 Renormalization of Mq 

At tree-level and in the DR Q renormalization scheme SU(2) l symmetry requires that the 
parameter Mq in the t and h mass matrices have the same value. This is, however, not the 
case at loop-level in the on-shell scheme due to different shifts <5M? in the t and in the b 



sectors |)2|, |93| . Therefore, for a combined treatement of stops and sbottoms, an appropriate 
renormalization procedure for M~ is necessary. 

Taking on-shell squark masses and mixing angles as input parameters, we define the on-shell 

M? as 

M%(q) = m 2 ^ cos 2 8g + m| 2 sin 2 9g — m| cos 2(3 (7| L — e q sin 2 8w) — m l , (3.37) 
see Eqs. (1.46) and (1.48). M|(t) and M|(t) are then related by 

M|(6) = M|(f) + 5M|(f) - 5M?(6) (3.38) 

with 

5M| (q) = bm\ cos 2 5 - + Sm 2 - 2 sin 2 0$ - (m 2 ^ - m? 2 ) sin 28 q 58 q - 2m q 8m q . (3.39) 

The underlying SU{2)l symmetry is reflected by the fact that the shift <5M~(t) — 5M^{b) is 
finite. 



4 DR = Minimal substraction with DRED. 
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3.4.3 Renormalization of m q A q 

The squark couplings to Higgs bosons involve the parameter A q . Therefore, also this parameter 
has to be renormalized at 0(a s ). We proceed like for the definition of the on-shell Mq and 
write the expression m q A q in terms of on-shell squark masses m qi and mixing angles 8 q , see 
Eqs. ([mi ) and ( |T30| ). We thus get [0 



5{m q A q ) = \ (<frn| - 5m 2 - 2 ) sm29 q + (m| - m\ 2 ) cos 29 q 59 q + /i {cot/5, tan/5} 5m q (3.40) 
where cot (3 (tan/5) has to be taken for q = t (b). 

The squark masses and couplings also depend on the parameters \x and tan/3. However, the on- 
shell fi and tan /5 are defined via electroweak processes. Hence /i and tan {3 are not renormalized 
at o{a s ). 
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3.5 Decays into Charginos and Neutralinos 

The tree-level amplitude for the decay q { — > q'xf is (see Eq. ( |1.113| ) and Fig. |3T4la): 

M°(q t - q'xf) = tgu(k 2 ) [k%P L + t;^). (3.41) 

The decay width at tree level is thus given by 

r°(ft - q'xf) = w 3 ([(4) 2 + (4) 2 ] * - ( 3 - 42 ) 

with X = m| -m^, — m? + . The 0(a s ) loop corrected decay amplitude is obtained by the shifts 

4 -4+«8 ,,) +*4 w +*4 (e) . ( 3 - 43 ) 
4^4+ 5 4 M + + ^ c) > ( 3 - 44 ) 

where the superscript v denotes vertex corrections, w wave function corrections, and c the shift 
from bare to on-shell couplings. The virtual corrections to the decay width, ST^ (a = v, w, c), 
can thus be written as 

g 2 K(m 2 -m 2 ,,m 2 + ) 

ffW «-^> - — s^r^ < 3 - 45 > 

[2 4 6^ + 2 4 54 (a) ) X - 4 m 3 ,m^ + (4 Sk^ + k% 6^) 
Analogously, one gets for squark decays into neutralinos 

2/2 2 2 \ 

a ' K(m~ ,m„,m~~ ) 

r°te-gx° fc ) = w- ( [(Qffc)2 + {b!k)2] X " 4a »^ m x°J . ( 3 - 46 ) 

*r<°>fe ,X2) = ( X m ' *' (3-47) 



X 



[2 a\ k 54^ + 2 bf k 5b% a) ) X - 4 m^o (<& <^ a) + ^ ) 



where X = m|. — m 2 —m 2 . In the following, we give our results for squark decays into charginos. 
The analogous expressions for decays into neutralinos are obtained by the replacements xf ~^ 
Xl q' - q, X -> X, 4 - at, *$ - &£, ^ - 5af fc (a) , and *4<°> - <«4 (a) . The shifts 



^4 j <^4 > etc. get contributions from gluon exchange, gluino exchange, and the four-squark 
interaction. As we will see, in the renormalization scheme used the contribution due to the 
four-squark interaction cancels. 
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Figure 3.4: Feynman diagrams for the 0(a s ) SUSY-QCD corrections to squark decays into 
charginos and neutralinos: (a) tree level, (b) and (c) vertex corrections, and (d) and (e) real 
gluon emission. For wave-function corrections see Figs. |3.1| and [372. 
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3.5.1 Vertex Corrections 



The gluonic vertex correction (Fig. [T4]b) yields 



SI 



5k 



a s 
3^ 



a s 
3^ 



(Ami + 2X)(C + d + C 2 ) + (2mi + + X)C X + B 



2m g im^+C 2 



k q 

ij f ' 



(4mJ + 2X)(C + d + C 2 ) + (2m|+ + X)d + B 



k q 



2m q >m^+C 2 



(3.48) 



(3.49) 



B , Co, dj and C 2 are the standard two- and three-point functions [fSSfl . In this case, -E> = 
£? (m-+, m? m^/) and C m = C m (m~., m^ + , m 2 ,,; A 2 , m~ , m 2 ,), where we follow the conventions 

of |86[. As usually, we introduce a gluon mass A for the regularization of the infrared divergence. 



The contribution to due to the graph of Fig. ^T4jc with a gluino and a squark q' n (n = 1, 2) in 
the loop is: 



m 3 - 



2a. 

3 7T 



'jn q ,a LR + m,a flL - m s a LL ) £ q nj + rn-+a RL k q n ^\Ci 



+ m q i 
+ m g 



{mq'OtRL - m q a LR + m g a RR )k q nj - m~+a LR £ q nj \ (d + d) 



(3.50) 



m q ,a RR - m q a LL + m g a RL )k q nj - m~+a LL £ q nj C + (XC 1 + B ) a RL k 



2 a, 



3 7T 



(m q >a RL + m q a LR - m- g a RR )k q nj + m~+a LR £ q nj d 



+ me 



{m q ia L R - m q a RL + m g a LL ) £ q nj - m-+a RL k q nj (d + d) 



(3.51) 



m q >a LL - m q a RR + m g a LR ) £ q nj - m~+a RR k q nj C + (XC X + B ) a LR 



91 . 

".I 



with 



aiL = (a LL ) in = Tllx Kii, a LR = (a LR ) in = K q x 1Z q n2 , 
a RL = (a RL ) in = Tl q 2 K q nl , a RR = (a RR ) in = TZ q 2 Tl q n2 . 

Here, B, = B„(m? t ,<,m^), and C m = t;^^^^ 



(3.52) 



3.5.2 Wave— Function Correction 



The wave-function correction is given by (i ^ i') 
ilSZf+SZufa)] £l j + 8Z Vi {q i )il j , 

: \5Zf + 6Z u (qj] k% + 5Z Vi {q>) k\ tj 



6k: 



(3.53) 
(3.54) 
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rL,R 



are the quark wave-function renormalization constants due to gluon and gluino exchange 
(Fig. |3.1|) . The squark wave-function renormalization constants Z n i{q^) stem from gluon, gluino, 
and squark loops (Fig. |3.2|). See Sect. for the explizit expressions. 



3.5.3 Renormalization of the Bare Couplings 

In order to make the shift from the bare to the on-shell couplings it is necessary to renormalize 
the quark mass as well as the squark mixing angle: 



= O q jn 5TZ q n + K% 50] 2 , 5kf ] = 0% SKI + 72* 50%, (3.55) 



5aT = A q kn 5TZ q n + Kh 5h q m , Sb%°> = B q kn 5K q n + n\ x 5h q Lk , (3.56) 
with 

sn , = f-mn0 9 cosU 
V — cost^ —smVq I 



50% = _ Vj2 5m h 50% = _ Uj2 5m b , (3.58) 

V 2 mw sin /3 y 2 mw cos (5 



and analogously for 5h q Lk and 5h q Rk according to Eqs. (1.120) - (1.123). 

For the renormalization of the squark mixing angle we use the scheme of [37] as explained in 



Sect. |3.4.1| . With this choice of 59g the squark contribution to the correction is zero: 5T^ W '^ + 



<5r(c,<?) = 0. Moreover, the off-diagonal contribution (i ^ j) of the gluino-quark loop in Fig. |3~2 
vanishes in this scheme. 



3.5.4 Real Gluon Emission 



The total virtual correction 5T virt = + + n*' is UV finite. In order to cancel also the 
infrared divergence due to A — > we include the emission of real (hard and soft) gluons, see 
Fig. [T3]d,e. STnatfe -> q'xj) = T(^ -> gq'xf) is: 



5^^ -> q'xj) 



g 2 g s 

67r 2 m^ 



(4) 2 + (4) 2 1 (A + J ) 

+ 2Z 



i.i • j ' - 1 - (3-59) 
to|.J o + w^Jn + (m| + ml, - m 2 -+)I 01 + I + I x 



where Z = (k q j) 2 + (ifj) 2 X — 4k q j£ q jm q >m^+. The phase space integrals J, I n , I nm , and I™ 

have (mg i ,m q ',m-,+) as arguments and are given in [BB]. An analogous expression holds for 

5r rea/ (g i -> gx°) = r(g i -> #<?Xfc)- The complete 0(a s ) correction 5T = 5T virt + 5T real is UV 
and IR finite. 
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3.5.5 Numerical Results 



Let us now turn to the numerical analysis. Masses and couplings of charginos and neutralinos 
depent on the parameters M, fi, and tan (3 (M' = | tan 2 6\yM). For the stop sector we use m^, 
mi 2 , cos 9f, [A, and tan f3 as input values. The sbottom masses and mixing angle are fixed by 
the assumptions = 1.12 Mg(t) and Aj, = A t . The other parameters are taken as explained 
in 



3.2.3 



Figure |3.5| shows the SUSY-QCD corrections 5T = T — T for the decay ti — > b Xi relative 
to the tree-level width as a function of for m~+ = 150 GeV, m^ 2 = 500 GeV, cos % = 0.65, 
and tan (3 = 3. In order to study the dependence on the nature of the chargino (gaugino- or 
higgsino-like) , we have chosen three sets of M and \x values: M <C (M = 163 GeV, /x = 500 
GeV), M > |/i| (M = 500 GeV, // = 163 GeV), and M ~ \p\ (M = p = 219 GeV). Near the 
threshold, the corrections strongly depend on m^; for larger mass differences this dependence 
is much weaker. We consider > 200 GeV: For a gaugino-like xf t ne correction is up to 
— 10% while for M ~ the correction typically amounts to —10% to —15%. The biggest effect 
is found for a higgsino-like chargino due to the large top Yukawa coupling; in this case we have 
ST/T° ~ —20% to —25%. The dependence on cos 6^ is shown in Fig. [J15] where we plot the 
0(a s ) corrected decay widths (full lines) together with the tree-level widths (dashed lines) for 
nif = 250 GeV and the other parameters as above. Again, one can see that the SUSY-QCD 
corrections can considerably change the decay width of t\. This is in particular true for the 
case \fi\ ^> M for which the correction is up to —40%. However, also for M ~ and M ^> 
one can have large corrections, especially if the tree-level width is very small. In our examples, 
the SUSY-QCD corrections to the ti bxt decay width are mostly negative. However, they 
can also enhance the tree-level width, see e.g. Figs. 2 and 3 of Ref. |f77| . 

As an example for stop decays into neutralinos we discuss the decay t<i — > t\\- Here we 
fix m^o = 100 GeV. Again, we take tan/3 = 3 and study three different scenarios: M 
(M = 208 GeV, // = 500 GeV), M > \y\ (M = 500 GeV, /i = 123 GeV), and M ~ (M = 
jj, = 230 GeV). Figure [T7] shows the SUSY-QCD corrections relative to the tree-level decay 
width as a function of m ( - , for m t - = 250 GeV and cos% = 0.65. In the shown range of m^ 2 , 
Sr/T° varies from 18% to 3% in case of a gaugino-like Xv m case °f a higgsino-like neutralino 
ST/T° = —26% to —37%, and in the mixed scenario 6T /T° = —46% to —26%. The dependence 



on cos #4 is shown in Fig. |3T8| . Analogously to Fig. [3]6] we here plot the 0(a s ) corrected decay 



widths (full lines) together with the tree-level widths (dashed lines) for m t - 2 = 500 GeV and 



the other parameters as in Fig. pO. Again, a striking effect can be seen for \fx\ M. The 



neutralino, chargino, and sbottom masses for the various (M, /i) values are listed in Tab. 3.1 . 



For completeness, we also briefly discuss sbottom decays into charginos/neutralinos for the 
parameters of Tab. |3.1| . In these scenarios, ~ 390 GeV, m~ b2 ~ 425 GeV, and cos^ ~ 1. 
Table |3.2| lists the SUSY-QCD corrected widths together with the relative corrections for all 
bi —* txj and hi — > bxt decays that are kinematically allowed. As can be seen, also for sbottom 
decays the SUSY-QCD corrections can be of either sign and amount to a few ten percent. 
Note, that in Tab. |3.2| tan/3 = 3. For larger values of tan/3 the corrections to the sbottom 
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Figure 3.5: SUSY-QCD corrections for the decay t\ bxf relative to the tree-level width as 
a function of m h , for = 150 GeV, m h = 500 GeV, cos0 t ~ = 0.65, M D = 1.12 Mq, A t = A b , 
and tan/3 = 3. Three scenarios are studied: (a) M = 163 GeV, /i = 500 GeV, (b) M = 500 
GeV, fj, = 163 GeV, and (c) M = fj, = 219 GeV. 



t — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




cos Qi 

Figure 3.6: Tree-level (dashed lines) and SUSY-QCD corrected (full lines) widths of the decay 
t\ —> bxf in GeV as a function of cos%, for m-+ = 150 GeV, = 250 GeV, mj 2 = 500 
GeV, Mf) = 1.12 Mq, A t = A b , and tan/3 = 3. Three scenarios are studied: (a) M = 163 GeV, 
fi = 500 GeV, (b) M = 500 GeV, fi = 163 GeV, and (c) M = fi = 219 GeV. 
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Figure 3.7: SUSY-QCD corrections for the decay t 2 — > txl relative to the tree-level width as 
a function of m t ~ 2 , for m^o = 100 GeV, m h = 250 GeV, cos0 t - = 0.65, M b = 1.12 Mq, A t = A b , 
and tan = 3. Three scenarios are studied: (a) M = 208 GeV, y. = 500 GeV, (b) M = 500 
GeV, fi = 123 GeV, and (c) M = // = 230 GeV. 
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Figure 3.8: Tree-level (dashed lines) and SUSY-QCD corrected (full lines) widths of the decay 
^2 — ► txl in GeV as a function of cos% for m^o = 100 GeV, = 250 GeV, m^ 2 = 500 GeV, 
Mjj = 1.12 Ma, A t = A b , and tan/3 = 3. Three scenarios are studied: (a) M = 208 GeV, 
y = 500 GeV, (b) M = 500 GeV, // = 123 GeV, and (c) M = \i = 230 GeV. 
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M 


n 




A. 2 


m x% 


m x\ 


m xt 


m xt 


m h 


m b 2 


163 


500 


78 


151 


503 


518 


150 


516 


384 


427 


500 


163 


136 


166 


262 


517 


150 


516 


391 


424 


219 


219 


94 


161 


223 


296 


150 


293 


388 


425 


208 


500 


100 


194 


502 


520 


193 


518 


385 


427 


500 


123 


100 


126 


259 


515 


112 


515 


391 


424 


230 


230 


100 


171 


234 


307 


161 


304 


388 


425 



Table 3.1: Neutralino, chargino, and sbottom masses for the scenarios discussed above (m^ = 
250 GeV, m h = 500 GeV, cos0 t ~ = 0.65, M b = 1.12 Mq, A t = A b , and tan (3 = 3). All values 
in [GeV]. 



decay widths become even more important since the bottom Yukawa coupling increases. 
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b 


i 


k 


(M, n) [GeV] 


channel 


T [GeV] 


sr/r° 


T [GeV] 


5T/T 


(163,500) 


A, 1 


0.11 


0.02 


0.2 


-0.02 




bx° 2 


1.01 


0.02 


0.09 


0.45 




/ v, 1 


0.83 


0.18 


0.09 


0.56 


(500,163) 




0.09 


-0.09 


0.01 


-0.54 




bxl 


(0.006) 


(-0.69) 


0.009 


-0.66 




bxl 


(0.003) 


(0.57) 


0.09 


0.02 




tXi 


2.03 


-0.11 


0.11 


-0.16 


(219,219) 


bxl 


0.22 


0.01 


0.16 


0.01 




bxl 


0.45 


0.05 


0.04 


0.05 




bxl 


(0.007) 


(-0.55) 


0.008 


-0.62 




bxl 


0.14 


0.21 


0.02 


-0.13 




txl 


0.76 


-0.01 


0.04 


-0.18 


(208,500) 


bx°i 

A. 1 


0.1 


0.05 


0.19 


-0.004 




bxi 


0.8 


0.06 


0.07 


0.58 




/V 1 


0.31 


0.68 


0.06 


0.76 


(500,123) 


bxl 

/V, 1 


0.09 


-0.14 


0.009 


-0.69 




bxl 


(0.009) 


(-0.66) 


0.01 


-0.66 




bxl 


(0.006) 


(0.42) 


0.09 


0.03 




tXi 


2.9 


-0.14 


0.11 


-0.16 


(230,230) 


bx\ 


0.2 


0.02 


0.16 


0.01 




bxl 


0.44 


0.06 


0.03 


0.08 




bxl 


(0.006) 


(-0.57) 


0.008 


-0.63 




bxl 


0.12 


0.25 


0.02 


-0.12 




tXi 


0.6 


0.02 


0.04 


-0.2 



Table 3.2: Sbottom decay widths at 0(a s ) and the relative corrections for the scenarios dis- 
cussed above {m h = 250 GeV, m h = 500 GeV, cos0 t ~ = 0.65, = 1.12 Mq, A t = A b , and 
tan/3 = 3). The values in brackets correspond to branching ratios below 1%. 
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3.6 Decays into W and Z Bosons 



At tree level the amplitude of a squark decay into a W ± or Z boson has the general form 



M°(q? - $V) = -ig c tjV (h + fc 2 )" e^kz), 



(3.60) 



with hi, k 2 , and k 3 the four-momenta of qf, q*j, and the vector boson V (V = W ± ,Z° 
respectively (Fig. |3]9|a). a and (3 are flavor indices. In the following we define rrii = 



1714 



ik 



TZ q jk for simplicity. Moreover, we shall use primes to explicitly 



distinguish between different flavors. With this notation the q^q^V couplings c^y are, see 
Eqs. ( [T95D and ( |Q6|) : 



1 



cos 8w 



7| L cos 2 6q — e q sin 2 9 W 



l\ L sin 2 9g — e q sin 2 9 



(3.61) 



COS 9g COS 6gt 
■ Sin 9q COS 9qi 



cos sin 9q> 
sin 6*,7 sin 6* 



'j 



(3.62) 



The tree-level decay width can thus be written as: 

g 2 (Q iy ) 2 /t 3 (m 2 ,m 2 ,my) 



167T my mf 



(3.63) 



The 0(a s ) loop corrected decay amplitude is obtained by the shift 

r . . T . _> r . ... _L A>) i X r ( w ) i A r ( c ) 



in fl3.60| ). The virtual correction to the decay width thus has the form 

r\ c\ q 2 K 3 (mf, m 2 , rrA) < , n \ 

J X7T mr \ my I ■> J 



(3.64) 



(3.65) 



The vertex correction stems from the five diagrams shown in Figs. |3.9| b-f. The gluon ex- 
change graphs of Figs. |3i|b-d yield 

r £? (m 2 , A 2 ,m 2 ) + B (m 2 , A 2 , m 2 



0C ijV 



o c ijV 
37T 



— 2 (mf + m 2 



rr%) (Co + Q + C 2 ) 



(3.66) 



with C m — C m (m 2 , my, m 2 ; A , m 2 , m 2 ). A gluon mass A is introduced to regularize the infrared 



divergence. The gluino-exchange contribution, Fig. 3.S e, gives 



8d 



(y,g) 

21Z 



a. 



2 m g niq (J| L — 2e q sin 2 9w) (Co + C\ + C 2 ) cos 2^ 



37T cos #w 

+ I\ L [ 2 m 2 Co + m 2 2 d + m 2 ^ C 2 + (m 2 - m 2 q ) (d + C 2 ) + £? 



(3.67) 



sin 20, 



9 I ' 
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Figure 3.9: Feynman diagrams for the 0(a s ) SUSY-QCD corrections to squark decays into 
vector bosons: (a) tree level, (b)-(f) vertex corrections, (g)-(i) real gluon emission. For wave- 
function corrections see Figs. |0] and 
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with C m = C m {va\ 2 ,m 2 z , m 2 qi ; m|, m 2 q , m 2 q ) and B = B (m z , m 2 , m 2 q ), for the decay q 2 -> q 1 Z°, 
and 



Sc. 



ijW 



V2 a, 

3 7T 



m~ a (C + Ci + C 2 ) (m g 7e i2 n' n + m q , K n K' j2 ) 

m? d + mj C 2 + m? (2C + d + C 2 ) + B 
- m q m q/ (Ci + C 2 ) 7£ i2 7^- 2 



(3.68) 



with C m = C m (m 2 , m 2 ,, m 2 ; m% m% w 2 ,) and B = B^m 2 ?, m 2 , m 2 ,), for the decay q { 
The squark loop of Fig. |315|f does not contribute because it is proportional to the four 
momentum of the vector boson. The total vertex correction is thus given by: 



Sr {v) 
0C ijV 



0C ijV + 0C ijV ■ 



(3.69) 



The wave— function correction is given by (i ^ ^ f) 



0C ijV 



5Zn(q?) + SZjjffi) c ijV + SZi>i(q°) c VjV + SZy^q?) dyy 



(3.70) 



with Z nm (q n ) the squark wave-function renormalization constants (see Sect. |3.3| ). 
The counterterms for the couplings are: 

^z = -^Loos29 q 59 q , 

and 



(3.71) 



Sr {c) -J- 

oc ijW -^2 



— sin 9q cos 6*g/ sin 9 q sin #g/ 

— COS 6q COS flg/ COS 6^ SU1 9g' 



+ 



— cos 9 q sin ^/ — cos 9 q cos 0g/ 



sin 6>g sin Q q i 



sin 6^ cos 9gi 



S9? 



(3.72) 



with 56^ as defined in Sect. 3.4.1 . For the decay q 2 — > <?i the counterterm fl3.71| ) completely 
cancels the off-diagonal wave function corrections (i ^ j in Fig. |3.1|) . In case of q { — > gj- 
the contribution of the squark bubble in Fig. |3.1| is cancelled. Thus, in both cases the total 
squark loop contribution to the correction is zero, <yr® = o. 



Finally, we add the emission of real gluons (Figs, ^JJg-i) in order to cancel the infrared 
divergence: 



5T real = r(gf - g$V) 

.2 „2 



3vr 2 



mi 



27- 



7 + h + m2 -^oo + I n + (m 2 + m 2 - m 2 ,) 7, 



'01 



(3.73) 
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with k = K(mj, m|, m v ); the phase space integrals /, I n , and I nm have (m^, rrij , my) as argu- 



3 

ments. 



Let us now turn to the numerical results. As the squark couplings to vector bosons 
depend only on the squark mixing angles, we just need the on-shell squark masses m^ 2 and 
mixing angles 9 g , and the gluino mass as input parameters. 



We first discuss the decay t 2 — > t\Z. Figure |3.10| shows the tree-level and the 0(a s ) SUSY- 



QCD corrected widths of this decay as a function of the lighter stop mass mi, for m t ~ 2 = 
650 GeV, cos 9 1 = —0.6, and m g = 500 GeV. SUSY-QCD corrections reduce the tree-level 
width by —11.7% to —6.8% in the range of m^ = 85 to 558 GeV. It is interesting to note that 
the gluonic correction decreases quickly with increasing while the correction due to gluino 
exchange varies only little with m^. In our example, 5T^ /T° = —4.5% and /T° = —7.2% 
at mi = 85 GeV; at = 550 GeV the gluonic correction is negligible whereas /T° is still 
-6.8%. 

Taking a closer look on the gluino mass dependence we find that the gluino decouples slowly. 
This is visualized in Fig. |3.11| where we plot the SUSY-QCD correction 8T = F — T of the 



decay t 2 — > t\Z relative to its tree-level width for m~ tx = 200 GeV, m^ = 650 GeV, and 
cos#£ = ±0.6. As can be seen, for large gluino masses 5T/T approaches ~ —3%. The negative 
spike at m g = 475 GeV is due to the t 2 — > tg threshold. Notice, that this threshold is less 
pronounced for cos^- > 0. 



The dependence on the stop mixing angle is shown in Figs. |3.12| and |3.13| . Figure |3.12| shows the 
tree-level width together with the 0(a s ) corrected width of t 2 — > t\Z as a function of cos#j, for 
mf a = 200 GeV, m^ 2 = 650 GeV, and = 500 GeV. With the tit 2 Z coupling proportional to 
sin 29 1 the decay width has maxima at cos% = ±^75 (maximal mixing) and vanishes in case of 
no mixing (cos^- = 0, ±1). SUSY-QCD corrections reduce the tree-level width by about —5% 



to —10%. The relative correction ST/T° for the parameters of Fig. |3.12| can be seen explicetly 



in Fig. 3.13| (black lines). In addition, we here also show the case = 1 TeV (gray lines). 



As the gluonic correction has the same 9% dependence as the tree-level width (in our example: 
Sr(9)/r° = -3%) the 9- t dependence in Fig. [TT| comes only from the correction due to gluino 
exchange. For m- g = 500 GeV and negative cos6^ the correction is of the order of —10% while 
for positive cos% it is roughly —5%. For m- g = 1 TeV the correction is about —3% to —4% 



with much less dependence on the stop mixing angle. Approaching cos^ = or ±1, <5r/r° 
diverges because the tree-level coupling c 2 \z vanishes while Sc 2 iz 7^ 0. In this case the decay 
width becomes of 0(a 2 s ). Note, however, that the appearance of this divergence, as well as the 
condition cos^- = {0, ±1}, is renormalization scheme dependent. 

The decay b 2 — > b±Z can be important for large tan/3 due to large mass splitting and mixing 
of bi,R- The SUSY-QCD corrections to this decay are similar to those to t 2 — > t\Z. However, 
the corrections due to gluino exchange are in general smaller and thus the dependece of 5T /T° 
on the sbottom mixing angle is very weak. 
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Figure 3.10: Tree-level (dashed line) and SUSY-QCD corrected (full line) widths of the decay 
£2 — ¥ t\Z as a function of m^, for m~ t2 = 650 GeV, cos0 4 ~ = —0.6, and m g = 500 GeV. 
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Figure 3.11: Tree-level (dashed line) and SUSY-QCD corrected (full line) widths of the decay 
I2 — > t\Z as a function of nig, for = 200 GeV, m ( - 2 = 650 GeV, and cos#f = ±0.6. 
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Figure 3.12: Tree-level (dashed line) and SUSY-QCD corrected (full line) widths of the decay 
t 2 — > t\Z as a function of cos% for mi = 200 GeV, m t ~ 2 = 650 GeV, and m g = 500 GeV. 




-1. -0.5 0. 0.5 

cos Qi 



Figure 3.13: SUSY-QCD correction relative to the tree-level width of the decay t 2 — > t\Z as 
a function of cos# 4 ~, for m~ tx = 200 GeV and m 4 ~ 2 = 650 GeV. The black (gray) lines are for 
m g = 500 (1000) GeV. 
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We next turn to the squark decays into W ± bosons. Here we discuss two special cases: 

(i) bi and b 2 decaying into a relatively light t\ plus W~ for small sbottom mixing (small tan f3 
scenario). In this case the mass difference of 61 and b 2 is expected to be rather small and thus 
the decays b 2 — > b\ (Z, h°, H°, A ) should be kinematically suppressed or even forbidden. 

(ii) A heavy t 2 decaying into a relatively light b\ plus W + for large sbottom mixing (large tan (3 
scenario) . 

Note, that here we have to take into account that (for a given value of tan/?) one of the 
parameters m^, m 4 ~ 2 , 9f, , , and d~ h is fixed by the others due to SU(2) L gauge symmetry, 



see Sect. 3.4.2. 



The decay widths of bi 2 — > t{W are shown in Fig. |3.14| as a function of the stop mixing 



angle, for m~ bi = 500 GeV, = 520 GeV, cos0 6 = -0.9 (at tree-level), m h = 200 GeV, 
nig = 520 GeV, and tan/3 = 3. The value of m% 2 is determined by the other parameters. Hence 
mi 2 varies from 533 GeV to 733 GeV depending on the stop mixing angle. Despite the larger 
phase space for the b 2 decay, the width of b 2 — ► t{W~ is smaller than that of 61 — > t{W~ 
because the W couples only to the "left" components of the squarks [b\ ~ b L and b 2 ~ b R for 
cos 6^ = —0.9). For | cos# f ;| > 0.1, the SUSY-QCD corrections change the tree-level widths by 
about —11% to +4%. For |cos6> t ~| — > 5T/T again diverges because the tree-level coupling 
vanishes. The corrections slowly decrease with increasing gluino mass e.g., 5T/T ~ —5% to 
— 1% for nig = 1 TeV and the other parameters as in Fig. |3.14j . The dependence of 5T /T° on 



the sbottom mixing angle is, in general, much weaker than that on the stop mixing angle (apart 
from a singularity in case of 6j = &r) . The overall dependence on the gluino mass is in general 
similar to that of the t 2 — > t 1 Z decay. However, the threshold effect at m g = m~ b . — m b is less 
pronounced. 



An example for large tan f3 is shown in Fig. |3.15| . Here we plot the tree-level and the SUSY 



QCD corrected widths of t 2 — > b±W + as a function of cos^- for m~ tx = 300 GeV, m 4 ~ 2 = 650 GeV, 
= 380 GeV, cos 6^ = —0.8, m g = 500 GeV, and tan/3 = 40. m^ 2 is calculated from the other 
parameters and thus varies from 615 GeV to 918 GeV. As expected, the decay width is maximal 
for t 2 = ti and vanishes for t 2 = t R . In the example chosen, the SUSY-QCD corrections are 
—2.4% to —4.7%. For m g = ITeV (and the other parameters as above), they are about —1% to 
— 1.5%. Again, there is almost no dependence on cos 0g apart from a singularity at cos6^ = 0, 
i.e. bx = b R . 
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Figure 3.14: Tree-level (dashed lines) and SUSY-QCD corrected (full lines) widths of the 
decays &i j2 — > t\W~ as a function of cos% for m~ hi = 500 GeV, m^ 2 = 520 GeV, cosO^ = —0.9, 
= 200 GeV, m g = 520 GeV, and tan (3 — 2. mi 2 is a function of the other parameters and 
varies with cos 9f. 
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Figure 3.15: Tree-level (dashed line) and SUSY-QCD corrected (full line) widths of the decay 
t 2 — > b 1 W + as a function of cos% for = 300 GeV, m t ~ 2 = 650 GeV, = 380 GeV, 
cos 9~ b = —0.8, rrig = 500 GeV, and tan/3 = 40. m~ b2 is a function of the other parameters and 
varies with cos 6*. 
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3.7 Decays into Higgs Bosons 



At tree level, amplitude and width of a squark decay into a Higgs boson, Fig. [3.16| a, are given 
by 



M°(q? -> f s H h ) = i{G?) v 



and 



\(G£)J 2 K(m 2 ,m 2 ,m 2 H ) 



167rm? 



(3.74) 



(3.75) 



where mi = and rrij = m-p. a and (3 are flavour indices and H k = {h ,H ,A ,H } 



fc = 1...4. For k — 1, 2, 3 we have of course a — j3 and i — 2, j — 1. For fc = 4 we have (gf , 



(ti, bj) or (bi,tj). The H^q^q" couplings (G^)^- are given in Sect. |1.6.4| . In the following, we 



will omit flavor indices when possible (flavor = a if not given otherwise). 
With the shift 



(G k h + 5{G k )^ + S(G k W + S(G k 



'ij 



the C(a s ) corrected decay amplitude can be expressed as 

/<m 2 ,m 2 ,mi ) 



5r (a) (gT _^ 



K{(G fc )*.5(G4 a) } (a = v,w,c). 



(3.76) 



(3.77) 



The vertex correction due to the gluon-squark-squark loop in Fig. 3.16 b is 



S(G k 



3tt 



G 



ijk 



B {m 2 , 0, m 2 ) + B (m 2 , 0, m 2 ) - B (m 2 Hk ,m 2 , m 2 ) + 2 X C 



(3.78) 



\(v,g) _ a s 
>ij 

with X = mf + m| — mfy and C = C (m 2 , m 2 lk ,m 2 ; A 2 , m 2 , m 2 ). Again, we introduce a gluon 
mass A to regularize the infrared divergence. 



The graph with the gluino-quark-quark loop in Fig. 3.16 c leads to 



r>Q C 

rrig cos zOq s< 

3 7T 



D / 2 2 2\ 



+ (4m 2 -m 2 )C 



(3.79) 



for the decays into h° and if (£ = 1, 2) 
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Figure 3.16: Feynman diagrams for the 0(a s ) SUSY-QCD corrections to squark decays into 
Higgs bosons: (a) tree level, (b)-(f) vertex corrections, (g)-(i) real gluon emission. For wave- 
function corrections see Figs. |3J] and 
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2a, 

3 7T 



| m q sin 26 q 



B, 



o(to? 2 ,to?,to 2 j ) 



*o(to|,to?,to; 



+ K -m?)q 



+ to 5 



B (m? 2 , to?, mj) + S (m? i; m|, m\) - m 2 A C ]} , (3.80) 



for the decay into A , and 



2 2^ 
m t ,m b ) 



(m^Au - m- g A 21 ) si + (m^v4 22 - m~ g A 12 ) sf B (m? , to?, mj«) 

J^l D A™2 _2 _2 



+ 
+ 
+ 



{m qf 3A 2 2 - m g A 21 ) + (m,^n - to 5 Ai 2 ) S4 J B o (m j ,m s ,m q0 ) 
m qf 3 (m 2 a - m- + to|) (A 22 S4 + An sf ) 
+ m g a (mjs - m 2 + to?) (A n s% + A 22 sf ) 
+ uig (mf + - TO^a - TOg/j) (A 2i S4 + A i2 sf ) 

— 2TOg TO ? a m g/ s (A.12 S4 + A.21 S4 ) C } (3.81) 



i,,„, - TZ° n 7lj m for the decay into a charged Higgs boson. 



with A n 

In Eqs. (|3.79|) to (|3.81|) C = Co(m 2 ,m 2 lk ,m 2 ;m 2 ,,m 2 a ,m 2 f3 ). The factors s q k are the Higgs 
couplings to quarks, see Eq. ( |1.100| ). 



The vertex correction due to the four-squark interaction in Fig. 3.16 d is 



8{G k tf - ° s ^ 1 



~7T~ ^2 Sjm G nmk B (TO^ , TTllp , TO 

d7r n,m=l,2 



with 



cos 26 q 



sin 26 q 



sin 26 q — cos 26 q 



(3.82) 



(3.83) 



The wave— function correction is given by 



6(G h 



j'j K^kJij' > 



i 7^ i' 
j * f 



(3.84) 



where the Z£ m are the squark wave-function renormalization constants for q a . 

We next fix the shift from the bare to the on— shell couplings 5(Gk)\f in Eq. (|3.76|) . From 
Eqs. ( |1.102| ) and ( |1.43j ) we get for the squark decays into h° or H° (£ — 1,2) 



S(G t ) 



«(c) 
21 



TV 8Gf (7^~) T + 571^ Gj (7^) T + TV Gf 5(11^ 



21 



COS 29 n 



5Gi 



21 



(G?)u ~ (GI) 



22 



69* 



(3.85) 
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with 8Gj obtained by varying Eqs. ( pL . 1Q5| ) and ( |1.108| ) e.g., 



5Gt 



5G b 2 



g ( Am t s a 5m t S(m t A t ) s a - /i c a 5m t 

2m w s P I 5(m t A t ) s a — fic a 5m t Am t s a 5m t 

g ( Am b c a 5m b 5{m b A b )c a - fis a 5m b \ 

2m w cp I 5(m b A b )c a — ns a 5m b Am b c a 8m b J 



(3.86) 
(3.87) 



with S/3 = sin/5, c@ = cos/?, s Q = sin a, c a = cos a, and a the Higgs mixing angle. SGf is 

PI by 

(3.88) 



obtained from Eqs. ( |3.86|) and (|3.87| ) by 

6G} = (5Gl with a -> a + ix/2) . 
For the couplings to the A boson we have explicitly 

*( G !)n = 1^7 [S(m q A q ) {cot/?,tan/?} + /^m 9 ], 
where cot (3 (tan (3) is for q = t (b). For the decay U — > bj H + {k = 4) we get 



6{G 



u (c) 
4)ij 



U l 8Gi {1Z 



b\T 



■lyiGj^sei-i-iyiG^se,, 



with 



5Ga 



2m b 5m b tan f3 + 2mt5mt cot j3 5(m b A b ) tan (3 + [i Sm b 
5(m t At) cot/5 + fi5m t 2(5m t m b + m t 5m b )/ sin 2/5 



(3.89) 



(3.90) 



(3.91) 



and analogously the expression for 6j — > tj H according to (|1.103|) . 

For the renormalization of m q A q and the squark mixing angels 9 q we use the prescriptions 



described in Sect. 3.4 



In order to cancel the infrared divergence we again include the emission of real gluons 
(Fig. |3~T6|e): 



5T real = T(qf -> H k g) 



a s \(G 



a\ 12 
k )ij\ 



?>TX 2 nrij 



I + h + m 2 I 00 + m 2 Iu + Xl 



01 



(3.92) 



Again, X = m 2 + m 2 — ra 2 Hk . The phase space integrals I n , and I nm have (m«, mj, rriH k ) as 
arguments. 

Let us now turn to the numerical analysis. As input parameters we use mt, m ( - 2 , cos^-, 
tan j3, fi, vriAi an d Trig. From these we calculate the values of the soft SUSY-breaking parameters 
Mg(t), Mq(1), M £)} and A tjb , taking M b = 1.12 Mq{1) and A b = A t for simplicity. 
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As a reference point we take m$ = 250 GeV, m ( - 2 = 600 GeV, cos% = 0.26 (9i ~ 75°), 
tan/? = 3, fi = 550 GeV, m A = 150 GeV, and m g = 600 GeV. This leads to m~ b = 563 
GeV Q, m- b2 = 627 GeV, cos0 6 = 0.99, A t)6 = -243 GeV, m h o = 100 GeV, m H o = 162 GeV, 
sina = —0.58, and m H + = 164 GeV. Thus t 2 can decay into t\ + (h°,H°,A°), and 5^2 can 
decay into t\H~ . 

We first discuss the parameter dependence of the widths of i 2 decays into neutral Higgs bosons 
by varying one of the input parameters of the reference point. Figure |3.17| shows the tree-level 
and the SUSY-QCD corrected widths of the decays t 2 — > t\ + (h°, H°, A ) as a function of m 4 - 2 . 
The relative corrections ST /T° = (T — r°)/r° are about —10% for the decay into h° and —9% 
to —62% for the decay into A . The corrections for t 2 — > tiH° are —9%, —45%, and +45% 
for mi 2 = 420, 670, and 900 GeV, respectively. The spikes in the corrected decay widths for 
mi 2 = 775 GeV are due to the t 2 — > tg threshold. The different shapes of the decay widths can 
be understood by the wide range of the parameters entering the Higgs couplings to stops. In 
the range m h = 300 GeV to 900 GeV, we have A t = 144 GeV to -889 GeV and sin a = -0.52 
to -0.73 (m h o = 81 GeV to 114 GeV, and m H o = 163 GeV to 170 GeV). 
Figure |3.18| shows the cos Qi dependence of the tree level and the corrected widths of t 2 — > 
ti + (h°, H°, A ) decays. Again the shapes of the decay widths reflect their dependence on the 
underlying SUSY parameters in a characteristic way. In particular we have A t = 1033, 183, 
-666 GeV and sin a = -0.748, -0.565, -0.726 for cos6 { = -0.7, 0, 0.7, respectively. Apart 
from the points where the tree-level decay amplitudes vanish the relative corrections range 
from -40% to 20%. 

In Fig. |3.19j we show the tree-level and the SUSY-QCD corrected decay widths as a function 



of m A . For m A = 100, 200, 300 GeV we have m h o = 85, 104, 105 GeV, m H o = 128, 207, 304 
GeV, and sin a = —0.87, —0.45, —0.37, respectively. The corrections to T°(t 2 — > tih°) range 
from -15% to -7% for m A = 100 GeV to 400 GeV. Those to T°(t 2 -> ^H ) are -50% to 
-22% for m A > 114 GeV, and those to T°(t 2 -> ^A ) are about -25%. 

As for the dependence on the gluino mass, the gluino decouples very slowly, as can be seen in 
Fig. In the range m~ g = 300 GeV to 1500 GeV <5r/r° varies from (-9%, -37%, -28%) to 



-7%, -16%, -14%) for the decays t 2 h + (h°,H°,A°), apart from the i 2 tg threshold 
at rrig = 425 GeV. 

As for the dependence on tan/5, we get T{t 2 -> ^/i ) = 2.68, 2.09, 1.42 GeV with <5r/r° ~ 
-10%, -7%, -5% for tan/3 = 3, 10, 30, respectively. Likewise, we get Y = 0.67, 1.61, 2.45 GeV 
with 5T/T ~ -27%, -19%, -17% for the decay into H° and T = 2.1, 2.64, 2.92 GeV with 
5T/T ~ -22%, -19%, -18% for the decay into A , respectively. 

Let us now turn to the sbottom decays. We start again from the reference point given above. 
For the decay 6i -> l x H- we get T = 3.88 GeV with 5T/T = -24%, and for the decay 
b 2 — > t\H~ we get V = 0.08 GeV with <5r/r° = +87%. As in our examples the width of 
the latter decay is usually quite small (because b 2 ~ 6r and t\ ~ t R ) we will discuss only the 
parameter dependence of the b\ decay. 



5 Notice that at tree level one has = 560 GeV because Mq = 558 GeV for both the t and b mass matrices. 



At 0(a s ), however, one gets Mg(t) = 558 GeV and M$(b) = 561 GeV. 
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Figure |3.21| shows the tree-level and the SUSY-QCD corrected widths of this decay as a func- 
tion of m^. The SUSY-QCD corrections are about —25%. Notice that at tree-level we have 
= 556 GeV to 566 GeV for = 85 GeV to 400 GeV, whereas at 0(a s ) we have 
m~ b = 561 GeV to 570 GeV. Therefore, the thresholds at tree-level and one-loop level are 
slightly different. 

The dependence on the stop mixing angle is shown in Fig. [3.22j for tan/3 = 3 and 10, and the 
other parameters as given above. (For | cos# t ~| > 0.72 the decay bi — > t\H~ is kinematically not 
allowed.) In case of tan/3 = 3, the SUSY-QCD corrections range from about —40% to 26%, 
with 5T/T > for cos^- < — 0.6. In case of tan/3 = 10 ST/T° is much larger. For cos#j > 0.5 
and tan/? = 10 we even get a negative corrected decay width. This is mainly due to a large 
contribution stemming from the term 5{m b A b ) ~ /itan/3 Suit of Eq. (|3.40Q and was already 
mentioned in [021. The same problem can occur for the decays b 2 ~~ > 61 + (h°,H°,A°) and 



t 2 — > biH + which may be important for large tan j3 (due to the large bottom Yukawa coupling 
and the large 61-62 mass splitting). 

We have also studied the dependence on itia- In the case tan/3 = 3 (10) we have found that 

<5r/r°(6i -> txH') 20% (-40%) for m A = 100 GeV to 285 GeV and the other parameters 

as given above. 

As for the dependence on the gluino mass, ST/T (bi — > t\H~) ranges from —26% to —14% 
(-47% to -39%) for m~ g = 300 GeV to 1500 GeV and tan/3 = 3 (10). 

In Ref. ]ST[] a numerical analysis was made for the decays t 2 — > t± + (h°, A ). Whereas we 
fairly agree with their Fig. 6 b for t 2 — > tiA° (apart from the fact that the spike is due to the 
t\ — > tg threshold), we find a difference of about 10% in both the tree-level and the corrected 
widths of t 2 — > t\h° (their Fig. 6 a). This may be due to a different treatment of the radiative 
corrections to the h° mass and mixing angle a [94|. For comparison, we show in Fig. p. 23 



our 



results on the t 2 —>■ t\h° decay for the parameters used in Fig. 6 a of Ref. In our notation 
these parameters are: My = 500 GeV, A t = fi = 300 GeV, tan/3 = 1.6, and uia = 400 GeV. 
Moreover, a s = 0.12. The shown range of is obtained by varying Mq from 200 to 430 
GeV. In this range, we get m^ 2 = 531 to 534 GeV, cos 9f ~ —1, m h o = 72 to 80 GeV, and 
sin a ~ —0.56. We have checked our results with two independent numerical programs. As a 
sidemark, we here also note that in the case of q. t — * q'jH Eq. (43) of Ref. [|8l]] is only correct 
if one takes the transposed couplings. 
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Figure 3.17: Tree-level (dashed lines) and 0(a s ) SUSY-QCD corrected (full lines) decay widths 
of t 2 -> ix + (h°, H°, A ) as a function of m ?2 , for m k = 250 GeV, cos0 t - = 0.26, [i = 550 GeV, 
tan/3 = 3, rriA = 150 GeV, and nig = 600 GeV. The grey area is excluded by the bound 
m h o < 90 GeV. 




Figure 3.18: Tree-level (dashed lines) and 0(a s ) SUSY-QCD corrected (full lines) decay widths 
of t 2 -> fx + (h°, H°, A ) as a function of cos 9 h for m h = 250 GeV, m h = 600 GeV, /i = 550 GeV, 
tan/3 = 3, ttia = 150 GeV, and m~ g = 600 GeV. The grey areas are excluded by the constraints 
given in Sect. |3.2.3| : Sp(t-b) > 0.0012 in (a), m h o < 90 GeV in (b), and unstable vacuum in (c). 



CHAPTER 3. SUSY-QCD CORRECTIONS TO SQUARK DECAYS 



92 




Figure 3.19: Tree-level (dashed lines) and 0(a s ) SUSY-QCD corrected (full lines) decay widths 
of t 2 -> h + {h°, H°, A ) as a function of m A , for m k = 250 GeV, m t ~ 2 = 600 GeV, cos 6- t = 0.26, 
[i = 550 GeV, tan/3 = 3, and m g = 600 GeV. The grey area is excluded by the bound 
m h o > 90 GeV. 
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Figure 3.20: 0(a s ) SUSY-QCD corrections (in %) to the widths of t 2 -> t x + H°, A ) as 
a function of m g , for = 250 GeV, m 4 ~ 2 = 600 GeV, cos^j = 0.26, \i = 550 GeV, tan/? = 3, 
and niA = 150 GeV. 
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Figure 3.21: Tree-level (dashed line) and 0{pt s ) SUSY-QCD corrected (full line) decay widths 



of b\ — > t\H as a function of mi, for mi 

50 GeV : 
one-loop level. 



600 GeV, cos0 t - = 0.26, fi = 550 GeV, tan/3 = 3, 



mA = 150 GeV, and m~ g = 600 GeV. The insert zooms on the different thresholds at tree- and 




Figure 3.22: Tree-level (dashed lines) and 0(a s ) SUSY-QCD corrected (full lines) decay widths 
of b\ — > i\H~ as a function of cos%, for = 250 GeV, m^ 2 = 600 GeV, /i = 550 GeV, 
niA = 150 GeV, m g = 600 GeV, and tan/3 = 3, 10. The grey area is excluded for tan/3 = 3 by 
the bound m^o > 90 GeV. 
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Figure 3.23: Tree-level (full line) and 0(a s ) SUSY-QCD corrected (dashed and dotted lines) 
decay widths of t% — > t\h° as a function of for the parameters used in Fig. 6 a of Ref. |81| . 
In our notation: M = 500 GeV, A t = fi = 300 GeV, tan/3 = 1.6, and m A = 400 GeV. The 
dashed line is for m~ g = 500 GeV and the dotted one for rrig = 1000 GeV. 
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3.8 Branching Ratios 



In this section we present the SUSY-QCD corrected branching ratios of stop and sbottom 
decays. We proceed as before, taking , mi 2 , cos % tan /3, fi, M, and as input parameters. 
The sbottom sector is fixed by M b = 1.12Mg(t) and A t 



i t — A h . Figure |3.24| shows the sbottom 
and Higgs boson masses obtained in this way as a function of the stop mixing angle for = 
250 GeV, m h = 600 GeV, fi = 500 GeV, tan/3 = 3, and m A = 200 GeV. 



For discussing the branching ratios of decays we take the parameter sets of Sect. [T5 
with m~+ = 150 GeV. Figures |3.25| a-c show the tree-level and the SUSY-QCD corrected 

branching ratios of t\ as functions of for = 600GeV, cos &i = —0.5, tan (3 = 3, and ttla = 
500 GeV. Moreover, M = 163 GeV, fi = 500 GeV in Fig. |3J5| a; M = 500 GeV, fi = 163 GeV 
in Fig. |3.25| b; and M = fi = 219 GeV in Fig. |3.25| c. If only the decay into b\t is kinematically 
allowed it has a branching ratio of practically 100%. Since the chargino decays further into 
X®qq' or ^{t^v the resulting signature is two acoplanar 6-jets plus additional jets and/or charged 
leptons plus large missing (transverse) energy, E™ ms > 2 rn^o. The decays ti — > txl (k = 1, ... 4) 
and ii — > b\t can a ls° have large branching ratios if they are kinematically allowed. Again, 
the decay cascades result in multi-jet/lepton signatures plus large missing energy. SUSY-QCD 



corrections change the tree-level branching ratios by O(±10%) in Figs. [3.25| a,b. However, they 
can also have larger effects. This can be seen in Fig. |3.25| c where the corrections go up to 
±30%. 



The dependence on the stop mixing angle is shown in Figs. |3.26| a-c. In order to have several 
decay modes open, we choose = 400 GeV and the other parameters as in Fig. |3.25| . The 
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Figure 3.24: Sbottom and Higgs boson masses as a function of cos^- for m$ = 250 GeV, 
m t - = 600 GeV, \i = 500 GeV, tan/3 = 3, and m A = 200 GeV; M b = 1.12Mo(t) and A t = A b . 
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branching ratios of t\ strongly depend on both the tjj-tn and the gaugino-higsino mixing. In 
Fig. |3.26j a, xt an d X12 are gaugino-like. For —0.5 < cos^ < 0.2 (Ji ~ t R ) the decay ti — > tx? 



dominates. Otherwise, t\ — > bxt is the dominant decay mode. BR(ti — > tx 2 ) has a similar 
cos^- dependence as BR^ — > bxt) but is phase-space suppressed. 

For higgsino-like xt an d X? 2 ) as m Fig- |3.26] b, the situation is quite different. In this case 
ti decays mainly into bxt unless t\ ~ t^. Moreover, BR^ — > tx®) > (<) BR^ — > tx 2 ) for 



cos flr < ( > ) -0.1. 



In the mixed gaugino-higgsino scenario of Fig. |3.26|c, ti decays mainly into bxt 2 with BR(t 



1 



bxt ) > 50% for cos^r > 0. The decays into neutralinos altogether do not exceed 30%. 

In Figs. |3.26 a and 3.26 b SUSY-QCD corrections change the tree-level branching ratios by 



C(±10%). In Fig. |3§c they go up to ±40%. 



Let us now turn to the branching ratios of t^- For the heavier stop mass eigenstate also 
decays into gauge and Higgs bosons may be important. We thus take m~ tl = 250 GeV, m t ~ 2 = 
600 GeV, rriA = 200 GeV, and tan (3 = 3 in this discussion. Figures |3.27| , |3.28| , and |3.29| show 



the cos6^~ dependence of the t 2 branching ratios for (M, /i) = (163, 500) GeV, (500, 163) GeV, 
and (219, 219) GeV, respectively. 

We first discuss the decays into charginos and neutralinos: In the case M <C BR(t 2 — > bxt) 
and BR(t 2 — > tx 2 ) are maximal for i 2 ~ ti (cos6* t ~ ~ 0). In the opposite case <C M, i 2 — Ir 
mainly decays into bxt an d h — ti mainly decays into txi 2 - For relatively small M ~ all 
decays into bxt 2 an d ^X? 4 are kinematically allowed. The branching ratios show an intricate 
dependence on the stop mixing angle. In particular, BR(t 2 — > tx 2 3 4) is maximal for t 2 ~ 
and BR(t 2 - bxt) > (<) BR(t 2 -> &*+) for cos% > (< ) 0.1. 

As for the bosonic decays, t 2 — > ii-Z is important in case of large t^-tn mixing. The decays 
into Higgs bosons are, quite generally, important for large values of \i and/or A t . The stop- 
Higgs couplings involve combinations /x ± A t . Since At = (m| — m~ )sin20f/(2mi) + //cot/3 
the branching ratios of t 2 — > + if , A ) decays show a distinct dependence on the sign 
of cos% Moreover, the tit 2 couplings to h° and if vanish in case of maximal mixing. In our 
analysis, the sbottom masses also vary with cos% For | cos# 4 ~| > 0.5 the decay t 2 — > b{W + is 
kinematically allowed. Its branching ratios go up to 40%. For | cos#r| > 0.8 also the decays 
into b 2 W + and biH + can be relevant. 



We have already noticed several times that the \i parameter plays a special role. Figures |3.30 



and [OT] show BR(t 2 ) a s a function of \i for m h = 250 GeV, m h = 600 GeV, M = 200 GeV, 



tan/3 = 3, and m A = 200 GeV. In Fig. cos0 t - = 0.4. For small |/i| ( < 400 GeV), decays 



into neutralinos dominate. For increasing |/i| decays into Higgs bosons become important 
or even dominant. Notice the dependence of BR(t 2 — > t\ + h°, H°) on the sign of fi. For 
300 GeV < ji < 800 GeV the decay i 2 — > t\Z has the largest branching ratio. The decays 
h —>■ bi + {W + , H + ) are kinematically not allowed. 

Figure pT30| has to be compared with Figure |3~3~T| where cos^ = 0.7. Here the decays into 



t\h° and t\H° have branching ratios below 1%; only the decay into t\A° can still have a large 
branching ratio. Moreover, the decay into b\W + is possible in this case and has a branching 
ratio of about 10% to 30%. The decay into t\Z is even more important than in Fig. |3. 30 . 
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Figure 3.25: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios of 
t± decays as a function of for m 4 ~ 2 = 600 GeV, cos % = —0.5, tan (3 = 3, htla = 500 GeV, (a) 
M = 163 GeV, /i = 500 GeV, (b) M = 500 GeV, \i = 163 GeV, and (c) M = /i = 219 GeV. 
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Figure 3.26: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios of 
ti decays as a function of cos % for = 400 GeV, mi 2 = 600 GeV, tan ft — 3, = 500 GeV, 
(a) M = 163 GeV, fi = 500 GeV, (b) M = 500 GeV, fi = 163 GeV, and (c) M = fi = 219 GeV. 
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For smaller m^-mx mass splitting and larger tua the bosonic decays are of course less important 
or kinematically forbidden. However, if is large the t 2 decays into t\Z, tih°, and biW + can 
still have large branching ratios. 

In all cases shown, SUSY-QCD corrections are important for precision measurements. 

Typical signals of ti decays into bosons and of those into fermions are shown in Table 1 of 
Ref. pa] . In principle, the final states of both types of decays can be identical. For example, 
the final state of the decay chain 

t 2 h + 0°, H°, A or Z) (txl) + (66) -> (bqq'xi) + (**) (3-93) 
has the same event topology as that of 

h^t + xIza - t + ((h°, H°, A or Z) + x°i) - t + (66"x?) ^ (bqq') + (66x?). (3.94) 
Likewise, 

£~ 2 -> 6 li2 + (H + or H/ + ) -> (bxl) + (qq) (3.95) 
has the same event topology as 

h -> 6 + -> 6 + ((^ + or iy + ) + x?) -> 6 + (qq'Xi)- (3.96) 

However, the decay structures and kinematics of the two modes (3.93) and (3.94) ((3.95) and 
(3.96)) are quite different from each other, since the Xi is emitted from t\ and X234 (^1,2 an d 
X12)) respectively This could result in significantly different event distributions (e.g., missing 
energy-momentum distribution) of the t% decays into gauge or Higgs bosons compared to the 
decays into fermions. Hence the possible dominance of the former decay modes could have an 
important impact on the search for i 2) and on the measurement of the MSSM parameters. The 
effects of the bosonic decays should therefore be included in the Monte Carlo studies of ti and 
62 decays. 
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Figure 3.27: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios of 
t 2 decays as a function of cos 9$ for = 250 GeV, m ( - 2 = 600 GeV, M = 163GeV, \i = 500 GeV, 
tan (3 = 3, and tua = 200 GeV; (a) shows the decays into fermions and (b) the decays into 
bosons. 



CHAPTER 3. SUSY-QCD CORRECTIONS TO SQUARK DECAYS 



101 




-1. -0.5 0. 0.5 1. 



cosOt 




-1. -0.5 0. 0.5 1. 

cos6^ 



Figure 3.28: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios of 
t 2 decays as a function of cos 6^- for = 250GeV, m~ t2 = 600GeV, M = 500GeV, \i = 163GeV, 
tan (3 = 3, and tua = 200 GeV; (a) shows the decays into fermions and (b) the decays into 
bosons. 
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Figure 3.29: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios 
of ti decays as a function of cos^ for = 250 GeV, m 4 ~ 2 = 600 GeV, M = \i = 219 GeV, 
tan (3 = 3, and = 200 GeV; (a) shows the decays into fermions and (b) the decays into 
bosons. 
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Figure 3.30: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios 
of t 2 decays as a function of \i for = 250 GeV, m^ 2 = 600 GeV, cos^- = 0.4, M = 200 GeV, 
tan/3 = 3, and uia = 200 GeV. The grey area is excluded by the bound m-+ > 95 GeV. 




Figure 3.31: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios of 
^2 decays as a function of /i for cos^- = 0.7 and the other parameters as in Fig. |3.30[ The grey 
area around fi = GeV is excluded by the bound m-+ > 95 GeV; in the gray area on the left 
the approximate condition for vacuum stability is violated. 
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Last but not least we discuss branching ratios of sbottom decays. Figure |3.32| shows 
the tree level and the SUSY-QCD corrected branching ratios of b± and h 2 decays as a function 
of the stop mixing angle for m h = 250 GeV, m~ t2 = 600 GeV, M = 230 GeV 0, /i = 500 GeV, 
tan/3 = 3, and = 200 GeV. The corresponding sbottom and Higgs boson masses are shown 
in Fig. Ipl Moreover, we have m^o = {111, 215, 502, 522} GeV, m~+ = {214, 520} GeV, and 



rrig = 645 GeV. For b\, which is mainly a bi, bosonic decays are important where they are 



kinematically allowed. To be concrete, in Fig. 3.32 b\ — > t\W is the dominant decay mode in 



case of large stop mixing because the W only couples to the "left" components of the squarks. 
For t\ ~ t R the decay b\ — > t\H~ is the dominant one because (i) b\t^W ± coupling vanishes 
and (ii) the coupling to H ± favours q^q'^ combinations. 

The situation is quite different for b 2 which is mainly a 6r. Here the decay into bxi dominates 
in most of the cos#j range. Only for 0.4 < cos% < 0.9 it has a branching ratio below 50%. In 
this range the decay into tiW~ dominates. 



The /i dependence of BR(&i j2 ) is shown in Fig. |3.33| for = 250 GeV, m^ 2 = 600 GeV, 
cos^j = 0.4, M = 200 GeV, tan/3 = 3, and uia = 200 GeV. Again, bosonic decays are 
very important for b\ while for b 2 decays into neutralinos are dominant in a wide range of \i. 
BR(6i,2 — * t\H~) increases with increasing \y\ because this parameter directly enters the tbH ± 
couplings. 



In both Figs. p.32| and |3.33| , SUSY-QCD corrections change the individual branching ratios by 
up to ±50%. 

As for the signatures of sbottom decays, analogous arguments apply as for the stop decays: 
If bip decays into b\\ the signature is two acoplanar 6-jets plus large missing energy. In 
all other cases one has multiple jets (and leptons) plus missing energy. Again, bosonic and 
fermionic decays can lead to the same final states. However, the various decay modes might 
be disentangled by their different event distributions due to the different decay structures and 
kinematics. 



6 We choose a somewhat higher value of M such that sbottom decays into gluinos are kinematically not 
allowed. 
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Figure 3.32: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios 
of b\ and 62 decays as a function of cos6^~ for = 250 GeV, m t ' 2 = 600 GeV, M = 230 GeV, 
fi = 500 GeV, tan = 3, and m A = 200 GeV. 
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Figure 3.33: Tree level (dashed lines) and SUSY-QCD corrected (full lines) branching ratios 
of b\ and 6 2 decays as a function of \i for = 250 GeV, = 600 GeV, cos 6^- = 0.4, 
M = 200 GeV, tan/3 = 3, and rriA = 200 GeV. The grey area is excluded by the bound 



m-+ > 95 GeV. 
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Open Questions 

In the recent years, extensive studies have been done for both stop and sbottom production 
and decays. However, some questions still remain open: 



The radiative corrections to squark prodution in have not yet been calculated. It 

can be expected that — similar to the e + e~ case — SUSY-QCD corrections and photon 
radiation are important. Indeed, a first analysis has shown that photon radiation can 
reduce the total production cross section significantly 

Since Yukawa couplings have a crucial influence on the t and b phenomenology, also 
Yukawa coupling corrections to the production and decay processes should be studied. 
So far, this has only been done for b { — > txj |$4 and e + e~ — > q^- [96]. In both cases it 



has turned out that these corrections are important. 

For SUSY-QCD corrections to squark interactions with Higgs bosons in the on-shell 
scheme, the parameter A q has to be renormalized. As mentioned in Sect. [3.7| , for large 
tan/? the term 5(mbAb) ~ /xtan /3 Srrib can become extremly large, leading to negative 
decay widths at 0(a s ). The inclusion of Yukawa corrections may help solve this problem. 
Another possible solution might be to take the bottom mass running. However, this is 
difficult to accomplish in a consistent way in the on-shell scheme. 
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